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ABSTRACT 


A  field  investigation  of  inner  shelf,  surf-zone,  and  beach 
processes  and  associated  morphodynamic  changes  was  carried  out  on  and 
seaward  of  Eastern  Beach,  on  the  Gippsland  Coast,  in  Eastern  Bass 
Strait,  Australia.  The  beach  is  part  of  a  long  straight  beach  system 
fronted  by  a  wide  continental  shelf.  Frequent  passage  of  gales  through 
Bass  Strait  results  in  relatively  high  energy  but  variable  and  rapidly 
changing  deepwater  wave  conditions.  Time  series  of  benthic  currents 
(u  and  v)  were  recorded  at  depths  of  20m  and  10m  over  the  nearshore 
profile  and  from  different  regions  of  the  surf  zone. 

Over  the  nearshore  profile,  wave-induced  oscillatory  flows 
dominate  to  depths  of  20m  and  are  responsible  for  agitating  and 
suspending  sediments.  Despite  the  microtidal  range,  reversing  shore- 
parallel  currents  and  superimposed  wind  driven  currents  are  second  in 
importance  and  determine  the  fate  of  sediments  suspended  by  waves. 
Frictional  dissipation  of  the  obliquely  incident  waves  over  the  rough, 
rippled  and  low-gradient  nearshore  profile  is  particularly  pronounced 
and,  as  a  result,  bottom  orbital  velocities  measured  at  10m  depth  were 
slightly  weaker  than  those  measured  at  20m.  Of  fundamental  importance 
to  surf  zone  processes,  nearshore  energy  dissipation  by  bed  friction 
causes  not  only  a  reduction  in  average  breaker  height  but  also  signifi¬ 
cantly  reduces  the  temporal  range  of  breaker  conditions. 

The  morphodynamic  state  of  the  surf  zone  and  beach  remains  inter¬ 
mediate  between  the  reflective  and  dissipative  extremes  and  is 
characterized  for  most  of  the  time  by  pronounced  longshore  bar-trough 
or  rhythmic  bar  and  beach  topography.  A  relatively  deep  (h=  2m)  trough 
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separates  a  dissipative  bar  from  a  steep  reflective  beach  face.  In 
addition  to  the  oscillations  at  incident  wave  frequency  which  are 
dominant  throughout  the  surf  zone,  important  roles  are  played  by 
surf  zone  standing  waves  at  subharmonic  and  infragravity  frequencies 
and  by  longshore  currents  and  rips.  Vertical  profiles  of  longshore 
current  velocity  within  the  trough  revealed  velocities  as  high  as 
0.25m  sec  1  to  within  0.075m  above  the  bed  and  were  not  at  all  loga¬ 
rithmic.  The  profiles  are  explained  in  terms  of  the  large  near-bottom 
addy  viscosity  caused  by  wave  induced  vortices  over  a  rippled  bed. 

Direct  measurement  of  vertical  profiles  of  suspended  sediment  concen¬ 
trations  showed  that,  despite  the  presence  of  longshore  currents, 
sediment  suspension  can  be  well  modelled  in  terms  of  entrainment  by 
waves  alone. 

Although  temporal  variability  of  deepwater  wave  conditions  is 
very  large,  the  substantial  reductions  in  breaker  variability  caused 
by  increased  nearshore  dissipation  of  larger  waves  also  reduces  the 
range  of  morphodynamic  variability  of  the  beach  and  surf  zone.  The 
longshore-bar-trough  and  rhythmic-bar-and-beach  states  represent  the 
modal  states  and  other  states  rarely  develop.  However,  profile  mobility 
is  considerable.  Most  of  this  mobility  is  expressed  by  onshore-offshore 
migrations  of  the  bar  and,  to  a  lesser  extent,  by  longshore  migrations 
of  rhythmic  features.  The  mobility  of  the  subaerial  beach  face  is 
comparatively  low  and  this  is  attributable  to  the  relative  constancy 
of  the  shore-break  amplitude.  The  Eastern  Beach  model  is  compared  to 
an  analog  at  Duck,  North  Carolina,  U.S.A. 
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INTRODUCTION 


It  is  apparent  that,  on  the  southeastern  coast  of  Australia  at 
least,  the  onshore-offshore  and  alongshore  exchange  of  sediments 
is  not  confined  to  the  surf  zone  but  extends  out  onto  the  inner 
continental  shelf  to  depths  of  over  20  metres.  In  addition, 
processes  operating  over  the  inner  continental  shelf  partially 
predetermine  the  nature,  intensity,  and  variability  of  the  processes 
which  operate  within  the  surf  zone.  In  turn,  the  complex  processes 
of  the  surf  zone  are  fundamental  in  effecting  the  transfer  of 
sediments  between  the  shelf  and  the  sui  aerial  beach.  The  problem 
of  beach  and  nearshore  morphodynamics  is  thus  one  which  involves 
the  entire  Coastal  Boundary  Layer  (Fig.l).  The  Coastal  Boundary 
•yer  is  defined  in  this  context  as  the  region  in  which  oceanic 
processes  are  measurably  modified  by  the  presence  of  the  coast. 

The  seaward  limit  is  approximately  the  50  metre  depth  contour  and 
the  landward  limit  is  the  beach.  In  this  report  we  are  concerned 
with  morphodynamic  aspects  of  the  Coastal  Boundary  Layer  including 
the  region  seaward  of  the  surf  zone,  the  surf  zone,  and  the  beach. 

All  deterministic  models  aimed  at  predicting  changes,  either 
erosional  or  accretional,  of  beaches,  nearshore  sea  beds,  or  any 
other  morphologic  feature  are  underlain  by  this  axiom:  morphodynamic 
changes  result  from  gradients  in  sediment  transport.  Sediment 
transport  depends  in  part  on  the  instantaneous  or  time  averaged 
shear  stresses,  X,  which  disturb  sediment  or  place  it  in  suspension 
and  in  part  on  the  ret  drift  pattern,  V,  which  advects  the  entrained 
sediment.  Fundamental  to  understanding  sediment  transporting 
processes  and  the; -  gradients  which  produce  either  accretion  or 
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erosion  is  an  appreciation  of  the  roles  played  by  different  modes 
of  fluid  motion  in  contributing  to  T  and  V  within  the  benthic 
regime  over  different  regions  of  the  Coastal  Boundary  Layer. 

Seaward  of  the  surf  zone,  we  must  understand  the  relative  roles 
played  by  interacting  waves,  tidal  currents,  and  other  currents 
in  determining  the  fate  of  nearshore  sediments.  In  the  presence 
of  complex  bar-trough  surf  zone  topography  the  process  signature 
is  even  more  complex  and  includes  the  effects  of  breaking  incident 
waves,  standing  waves  and  edge  waves,  and  wave-generated  longshore 
currents  and  rips. 

Our  field  studies  of  Australian  beaches  and  surf  zones  over  the 
past  few  years  have  shown  that  nearshore  and  surf  zone  process 
signatures,  associated  sediment  transport,  and  beach  response  are 
all  mutually  related  and  are  strongly  dependent  on  coastal  environ¬ 
mental  conditions  as  well  as  on  the  predominant  modal  state  of  the 
beach  and  surf  zone  (Wright  et  al  1979  a  and  b;  Short  1979  a  and  b; 
Wright  1981,  1982;  Wright  et  al  1981,  1982;  Wright  et  al  1982; 
Wright  and  Short  in  press) .  Depending  on  local  environmental 
conditions,  sediments,  inner  shelf  processes  and  antecedent  wave 
conditions,  beaches  and  surf  zones  may  be  dissipative,  reflective, 
or  in  one  of  four  intermediate  states.  The  intermediate  states, 
particularly  states  characterized  by  pronounced  longshore  bar  and 
trough  and  crescentic  bar  topographies  with  rip  circulation  cells 
are  typically  the  most  dynamic  in  terms  of  morphologic  mobility 
and  exhibit  the  most  complex  sets  of  processes  (Short,  1981; 

Wright  1981;  Wright  and  Short  in  press) . 

The  Gippsland  coast  of  Victoria  bordering  Bass  Strait  in  the 


extreme  southeastern  corner  of  Australia,  exhibits  an  uninterrupted 
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sandy  beach  (Ninety  Mile  Beach) ,  145  kilometres  in  length  fronted 
by  a  comparatively  wide,  petroleum  rich,  continenta  ^helf.  Ninety 
Mile  Beach  and  its  associated  nearshore  system  face  oughly  toward 
the  south  and  into  Bass  Strait  (Fig. 2)  through  which  the  frequent 
passage  of  mid-latitude  storms  and  gales  produce  rapidly  changing 
deep  water  wave  conditions.  This  highly  variable  and  relatively 
high  energy  wave  climate  combines  with  tidal  and  wind  driven 
currents  to  produce  a  complex  set  of  sediment-transporting  processes 
but  is  also  filtered  appreciably  by  the  wide  nearshore  profile  to 
reduce  the  temporal  variability  in  beach  morphology.  As  a 
consequence  of  the  regional  environmental  conditions  and  sediment 
properties,  the  modal  morphodynamic  state  of  the  beach  and  surf 
zone  tends  to  remain  intermediate  between  the  dissipative  and 
reflective  extremes  and  is  characterized  by  pronounced  bar-trough 
and  crescentic  bar  topography.  The  region  thus  offers  an  excellent 
laboratory  for  studying  the  processes  which  operate  in  the  presence 
of  this  globally- important  morphodynamic  state  as  well  as  the 
processes  which  operate  seaward  of  the  surf  zone  and  predetermine 
the  probability  that  this  surf  zone  state  will  exist. 

The  Gippsland  coast  provides  us  with  a  case  study  for 
illustrating  some  of  the  processes  which  operate  within  and  sea¬ 
ward  of  intermediate  bar-trough  and  rhythmic  bar  and  beach  type 
surf  zones.  The  functions  of  the  report  are  twofold:  (1)  to 
present  the  results  of  a  set  of  experiments  carried  out  on  and 
seaward  of  Eastern  Beach  near  the  eastern  end  of  Ninety  Mile  Beach 
in  eastern  Bass  Strait,  and  (2)  to  use  this  "case  study"  as  a 
medium  for  illustrating  the  application  of  recently  developed 
ideas  and  discussing  some  processes  of  general  importance.  To 


a  large  degree,  the  emphasis  is  on  principles  rather  than  place 
We  will  discuss  the  relative  magnitudes  of  several  important 
processes  at  different  positions  over  the  nearshore  zone  and  in 
the  surf  zone;  draw  some  inferences  about  nearshore  sediment 
transport;  and  discuss  the  beach  variability  and  mobility  that 
result  from  these  processes. 


THE  GIPPSLAND  COASTAL  ENVIRONMENT 


Coastal  Configuration  and  Geology 

The  Gippsland  coast  extends  440  km  along  an  arc  from  Wilsons 
Promontory  to  Cape  Howe  at  the  New  South  Wales  border  (Fig. 3). 

North  from  Snake  Island  the  coast  consists  of  an  almost  continuous 
364  km  sand  beach-barrier  system  facing  Bass  Strait.  The  low 
energy,  tidal  dominated  Shallow  Inlet  region  occupies  the  south¬ 
western  corner.  North  of  the  first  34  km,  exposure  to  wave  energy 
has  produced  a  continuous  beach,  Ninety  Mile  Beach,  extending 
from  Shallow  Inlet  to  Red  Bluff  a  distance  of  144  km  or  90  miles. 
Beyond  Red  Bluff  the  beach-barrier  continues  almost  unbroken  for 
another  132  km  to  Ram  Head.  Only  the  mouth  of  the  Snowy  River  and 
Cape  Conran,  Pearl  Point,  Clinton  Rocks  and  Point  Hicks  interrupt 
the  curve  of  the  coast.  Beyond  Ram  Head  the  coast  abruptly  trends 
northeast  resulting  in  initially  low  wave  energy  and  a  sediment 
deficient  coast.  Rocky  coastline  dominates  17  of  the  next  32  km 
with  sand  increasing  in  abundance  toward  Gabo  Island.  At  the  New 
South  Wales  border  a  wide  sand  plain  and  transgressive  coastal 
dunes  dominate  the  morphology.  In  total,  92%  of  the  336  km  between 
Shallow  Inlet  and  Cape  Howe  is  composed  of  sand  beaches  backed  by 
barriers  and  dunes.  Changes  in  coastal  orientation  along  the  arc 
(Fig. 3)  result  in  systematic  variation  in  levels  of  wave  energy 
and  onshore  wind  velocities,  which  in  turn  control  the  beach  and 
dune  sediment  supply  and  ultimately  the  form  of  the  beach-barrier 
system.  In  addition,  the  overall  orientation  relative  to  the 
dominant  west  and  southwest  waves  and  winds  produces  a  pronounced 
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northeastward  increase  in  wave  energy  and  associated  littoral 


drift. 

The  coast  lies  in  a  Tertiary  depositional  basin  called  the 
Gippsland  Basin  (see  Franklin  and  Clifton,  1971  for  a  description 
of  the  stratigraphy) .  It  is  bordered  on  the  north  and  west  by 
Palaeozoic  rocks  of  the  Eastern  Highlands,  and  in  the  southwest 
by  Cretaceous  rocks  of  the  Gippsland  Highlands  and  the  high  granite 
of  Wilsons  Promontory.  The  eastern  boundary  of  the  basin  lies  at 
the  edge  of  the  continental  shelf.  Because  of  the  present  high 
sea  level  the  shoreline  lies  approximately  in  the  middle  of  the 
basin.  East  of  Red  Bluff  the  east-west  trend  of  the  Eastern 
Highlands  dominates  the  coast  and  controls  the  coastal  orientation. 
South  of  Red  Bluff  the  present  coastline  lies  along  the  axis  of 
basin  tilting  with  slight  uplift  landward  and  extensive  downwarping 
seaward  across  the  shelf  with  sediments  exceeding  2000  m  in 
thickness.  The  Basin  is  composed  of  Tertiary  fluvial  and  marine 
deposits.  Uplifted  marine  deposits  outcrop  at  numerous  locations 
across  the  subaerial  plain.  During  the  Quaternary  at  least  three 
major  barrier  systems  were  built:  the  Prior,  Inner  and  Outer 
Barriers  (Bird  1978) ,  the  Outer  representing  the  Holocene  barrier 
system. 

Evidence  suggests  that  periodic  high  wave  energy  and  pronounced 
north-eastward  littoral  drift  has  delivered  Holocene  beach-barrier 
sediments  both  from  offshore  where  they  have  eroded  down  to 
Tertiary  sands  and  gravels,  and  from  alongshore,  resulting  in 
erosion  of  Pleistocene  deposits  between  Seaspray  and  Reeves  Beach. 
The  littoral  drift  has  resulted  in  continual  Holocene  inlet  and 
river  mouth  migration  northward  along  the  coast.  On  Ninety  Mile 
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Beach  numerous  abandoned  recurved  inlet  spits  and  flood  tide  deltas 
extend  from  Stockyard  Hill  for  16  km  to  Bunga  Arm,  and  Bunga  Arm 
itself  is  the  result  of  eastward  spit  growth  leaving  an  abandoned 
Holocene  barrier  abutting  Boole  Poole  Peninsula.  Similar  eastward 
spit  growth  and  truncation  was  a  characteristic  of  Cunningham  Arm 
before  the  construction  of  training  walls  at  Lakes  Entrance  (Bird, 
1978) . 


Climate  and  Wind  Regime 

The  Gippsland  coast  lies  between  latitudes  37°  34 'S  and  39°  08' S 
which  places  it  near  the  southern  margins  of  Australia's  dominant 
pressure  system,  the  subtropical  anticyclonic  belt,  and  in  the  path 
of  the  mid-latitude  cyclones.  Superimposed  on  the  general  westerly 
airflow,  stronger  westerly  and  southwesterly  winds  (locally  called 
gales)  accompany  the  passage  of  fronts.  Although  these  fronts 
occur  throughout  the  year,  they  are  most  frequent  and  their  effects 
most  intense  during  the  winter  months.  Following  passage  of  fronts, 
the  winds  usually  moderate  and  turn  south  or  southwest.  Northeast 
and  easterly  winds,  which  are  subdominant,  prevail  during  interven¬ 
ing  periods  of  high  pressure.  The  winter  period  of  most  frequent 
frontal  passage  is  also  the  period  of  greatest  rainfall.  Figure  4 
shows  the  mean  annual  distribution  of  temperature  and  rainfall  at 
the  two  ends  of  the  region.  Wind  roses  for  Tidal  River,  Lakes 
Entrance,  and  Gabo  Island  are  shown  in  Figure  5.  Additonal  wind 
data  can  be  found  in  Appendix  A. 

Winds  from  the  westerly  sector  (SW,  W,  NW)  dominate  at  ail 
stations.  However,  there  is  considerable  variation  along  the  coast, 


presumably  reflecting  local  topographic  effects.  There  is  also  a 
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Figure  5.  Wind  Roses  for  Tidal  River,  Lakes  Entrance  and  Gatoo  Island 
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significant  gradient  in  wind  velocity  between  Tidal  River  and 
Gabo  Island.  At  Tidal  River,  which  is  partially  sheltered  the 
"modal"  wind  speed  is  between  1-5  km  hr  ^  and  only  12%  of  the 
winds  exceed  21  km  hr  1,  at  the  more  exposed  locality  of  Lakes 
Entrance,  the  "modal"  wind  speed  is  about  10  km  hr  1  and  speeds 
of  21  km  hr  1  are  exceeded  for  20%  of  the  time.  The  strongest 
winds  prevail  at  Gabo  Island  which  is  fully  exposed:  there  the 
modal  speed  is  20  km  hr  1  which  is  exceeded  for  50%  of  the  time. 
At  all  three  locations,  the  strongest  winds  are  associated  with 
westerly  gales.  The  data  in  Appendix  A  suggests  that  diurnal 
seabreezes  also  contribute  significantly  to  the  wind  regime. 

At  all  stations  there  is  a  significant  increase  in  onshore 
velocities  between  0900  and  1500  hours. 

Wave  Climate 

The  coastal  wave  regime  of  the  region  is  controlled  by  four 
factors.  The  first  is  that  the  prevailing  westerlies  blow 
generally  offshore  to  the  south  of  Lakes  Entrance  and  alongshore 
to  the  east  whereas  the  southwesterly  winds  produce  obliquely- 
incident  waves.  Second,  fetch  lengths  from  the  west  and  south 
increase  eastward.  Third,  Wilsons  Promontory  and  offshore 
islands  north  of  Tasmania  partially,  but  not  wholly  reduce  the 
westerly  and  southwesterly  swell  reaching  eastern  Bass  Strait 
from  generation  centres  to  the  west  and  southwest  in  the  Southern 
Ocean.  Finally,  the  wide,  low-gradient  continental  shelf 
fronting  the  coast  reduces  the  heights  of  the  larger  waves  before 
they  reach  the  surf  zone.  This  effect  is  discussed  further  in 
the  next  section.  Overall,  the  wave  climate  is  thus  strongly 
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influenced  by  locally  generated  wind  waves  related  to  the  passage 
of  cyclones  and  gales.  Accordingly,  deep  water  wave  conditions  are 
subject  to  appreciable  and  rapid  temporal  variations. 

Although  deepwater  wave  conditions  have  been  recorded  from 
offshore  petroleum  platforms  for  several  years,  those  data  are 
proprietary  and  not  publicly  available.  We  have  used  some  of 
these  data  from  the  period  of  our  experiments  (with  permission  from 
Esso  Australia)  elsewhere  in  the  report;  however  we  do  not  have 
access  to  records  of  sufficient  length  to  describe  the  long-term 
wave  climate.  Hence,  we  base  our  description  here  on  two  sources: 

(1)  hindcasted  wave  statistics  for  western  Bass  Strait  for  1951- 
1965  (Spillane  et  al  1972);  and  (2)  visual  observation  data  from 
the  Gabo  Island  lighthouse.  Summaries  of  these  data  can  be  found 
in  Appendix  A.  Unfortunately,  both  sets  of  data  must  be  considered 
suspect,  at  least  quantitatively. 

Figure  6a  shows  the  average  annual  distribution,  by  direction, 
of  energy  flux  as  computed  from  the  hindcasted  wave  statistics  for 
western  Bass  Strait  (Spillane  1972) .  Figure  6b  shows  the  percentage 
of  occurrence  of  waves  arriving  from  different  directions  at  Gabo 
Island.  The  western  Bass  Strait  data  suggests  that  the  dominant 
energy  flux  arrives  from  the  southwest.  However,  the  Gabo  Island 
data  indicate  that  southerly  waves  dominate,  at  least  in  terms  of 
frequency  of  occurrence. 

Figure  7  indicates  the  probability  of  exceedance  of  wave  heights 
based  on  data  from  both  sources.  It  is  worthy  of  note  that  the 
hindcasted  data  from  western  Bass  Strait  and  the  sea  state  data  from 
Gabo  Island  show  reasonable  agreement.  Swell  height  data  from 
Gabo  Island  (which  are  suspect  owing  to  the  difficulty  of  obtaining 
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Figure  6.  Deepwater  Wave  Climate  of  Bass  Strait:  (a)  directional 
distribution  of  energy  flux  (watts  per  metre)  based  on 
hindcast  data  for  western  Bass  Strait;  (b)  frequency  of 
occurrence  of  waves  from  different  directions  based  on 
lighthouse  observations  at  Gabo  Island. 


PERCENT  EXCEEDANCE 


Figure  7.  Bass  Strait  wave  climate:  average  annual  exceedance  of 
deepwater  wave  heights,  H^.  Based  on  hindcasted  data 
for  western  Bass  Strait  and  lighthouse  observation  from 
Gabo  Island. 
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accurate  visual  estimates  of  swell  height  when  seas  are  superimposed) 
reveal  an  unusually  high  frequency  of  occurrence  of  high  swell.  If 
we  take  the  line  on  Figure  7  as  representative,  we  can  make  the 
following  generalizations  about  the  wave  climates  deepwater  waves 
exceed  1  metre  for  40%  of  the  time;  they  exceed  2  metres  for  20% 
of  the  time;  they  exceed  3  metres  for  4%  of  the  time;  and  they 
exceed  5  metres  for  nearly  2%  of  the  time.  From  the  western  Bass 
Strait  data  (Spillane  et  al  1972) ,  we  estimate  the  modal  wave 
(defined  as  the  set  of  conditions  for  which  the  product  of  energy 
flux  and  frequency  of  occurrence  is  at  a  maximum)  to  be  a  3  metre, 
8.5  second  southwesterly  wave.  The  "average"  (root-mean-square) 
wave  conditions  are  1.6  m,  8.5  seconds  from  the  data  of  Spillane 
et  al  (1972)  and  roughly  1.5  m,  10  seconds  from  the  Gabo  Island 
lighthouse  data.  The  "storminess"  of  the  Bass  Strait  environment 
is  reflected  in  the  hindcasting  statistics  which  indicate  that  the 
average  interval  separating  the  occurrence  of  waves  of  2  metres  or 
more  in  height  is  a  little  less  than  3  days. 

Tides 

The  tides  in  eastern  Bass  Strait  are  semidiurnal  with  a  moderate 
diurnal  inequality.  The  M2  tide  is  strongly  dominant  and  the  S2 
constituent  is  low  in  amplitude  so  that  neap-to-spring  variations 
are  modest.  The  Australian  Tide  Tables  indicate  mean  spring  tide 
ranges  of  1.0  m  at  Gabo  Island  and  2.1  m  in  the  Kent  Island  Group 
southeast  of  Wilsons  Promontory.  The  approximate  mean  spring  range 
at  Lakes  Entrance  is  only  0.9  m;  however,  this  estimate  is  based 
on  a  gaging  station  located  inside  the  entrance  and  is  thus  likely 
to  be  an  underestimate  of  the  "open  coast"  tides. 


METHODS 


This  report  is  based  largely  on  a  set  of  experiments  carried 
out  at  Eastern  Beach,  near  Lakes  Entrance  (Fig. 2)  over  a  relatively 
short  period  from  4  May  through  16  May,  1981,  supplemented  by  an 
additional  set  of  surveys  made  in  December  1981.  Despite  the 
limited  duration,  the  main  experiments  were  conducted  during  one 
of  the  stormiest  months  of  the  year  and  included  measurements 
obtained  during  two  gales  as  well  as  during  intervening  calm 
periods.  The  data  should  thus  be  reasonably  representative  of  the 
modal  conditions. 

Surveys  of  the  beach  and  surf  zone  were  carried  out  using  a 
conventional  theodolite  and  a  surveying  level.  A  specially  designed 
staff  permitted  the  staff  handler  to  operate  in  the  surf  zone  to 
water  depths  of  about  4  metres.  Surveys  were  extended  seaward  of 
the  surf  zone  using  a  Raytheon  Model  DE719  surveying  echo  sounder 
on  a  7  metre  workboat.  Horizontal  positions  for  bathymetric 
surveys  were  determined  by  theodolite  angles  from  shore.  At  the 
main  Eastern  Beach  experiment  site  a  survey  grid  was  established 
consisting  of  six  profiles  spaced  at  intervals  of  50  metres 
alongshore.  These  lines  were  surveyed  daily  over  the  experiment 
period.  Nearshore  sediment  samples  were  obtained  by  divers  and 
by  means  of  a  dredge  from  the  survey  boat.  Ripple  dimensions  on 
the  nearshore  bed  were  measured  directly  by  divers. 

Benthic  currents  seaward  of  the  surf  zone  were  recorded  from 
depths  of  10  metres  anc!  20  metres  using  a  pair  of  InterOcean  Model 
195  m  recording  electromagnetic  current  meters.  These  instruments 


employ  Marsh  McBirney  two-component  electromagnetic  current  sensors 
with  a  range  of  0- 3  m  sec  ^  and  a  precision  of  ±  2  cm  sec  \  They 
may  be  programmed  for  different  sampling  intervals,  burst  durations 
and  intervals  between  bursts.  In  this  case,  we  used  a  sampling 
interval  (At)  of  2  seconds,  a  burst  duration  of  15  minutes,  and  an 
interval  of  1  hour  between  the  start  of  each  burst.  The  instruments 
were  mounted  on  tripods  weighing  220  kilogrammes  with  sensors 
situated  at  elevations  of  1  metre  above  the  bed.  Figure  8  shows 
one  of  these  instruments  mounted  on  its  tripod.  Data  were 
recorded  on  magnetic  tape  cassettes  and  subsequently  read  into 
files  for  analysis  on  the  University  of  Sydney’s  Cyber  computer. 

Within  and  near  the  surf  zone,  strain-gage- type  pressure 
transducers  with  0-5  bar  range  absolute  were  used  to  measure 
pressure  which  was  converted  to  water  surface  elevation,  n. 
Horizontal  surf  zone  currents  were  measured  by  means  of  small,  low 
inertia  bidirectional  ducted- impel lor  flow  meters  (Bradshaw  et  al 
1978;  Sonu  et  al  1974) .  Flow  meters  were  mounted  orthogonally  to 
obtain  time  series  of  u  and  v.  Pressure  transducers  and  flow  meters 
were  mounted  together  on  masts  with  50  kilogramme  flat  bases  in 
various  arrays,  an  example  of  which  is  shown  in  Figure  9.  Several 
such  installations  were  deployed  simultaneously.  Signals  from  the 
sensors  were  transmitted  by  cables  to  a  mobile  van  housing  a  chart 
recorder,  multiplexing  digital  tape  recording  facilities  and  mini¬ 
computers.  The  field  monitoring  system  has  been  described  by 
Bradshaw  et  al  (1978) .  Sensors  were  sampled  simultaneously  at 
intervals  (At)  of  1-1.5  seconds.  Data  were  corrected  for  frequency 
response  characteristics  (Nielsen  and  Cowell,  1981) . 

Data  tapes  from  nearshore  and  surf  zone  instruments  were 
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t  r.msf erred  into  files  for  analysis  on  the  University  of  Sydney's 
Cyber  computing  system.  Analyses  included  computing  instantaneous 
and  t ime-average  vector  resultants  and  performing  spectral  and 
cross-spectral  analyses  to  determine  the  magnitudes  of  currents 
and  surface  oscillations  corresponding  to  different  frequencies 
and  forcings.  A  fast  Fourier  transform  was  used  to  compute 
spectral  estimates.  A  phase  relationship  of  90°  (t/2)  between 
water  surface  elevation,  n,  and  shore-normal  current,  u,  at  any 
given  frequency  was  used  as  a  test  for  the  existence  of  surf  zone 
standinq  waves  at  that  frequency.  Current  and  surface  oscillation 
amplitudes  related  to  different  frequency  bands  were  determined 
from  the  total  variance  within  that  band. 

Vertical  profiles  of  suspended  sediment  concentration  in  the 
surf zone  were  measured  using  the  suction  sediment  sampler  developed 
by  Nielsen  (Nielsen  and  Green  in  press).  This  device,  shown  in 
Figure  10  has  been  extensively  field  tested  and  proven  to  be  the 
most  effective  system  for  obtaining  concentration  profiles 
presently  in  use.  Water  and  sediment  enters  by  way  of  thin 
me'al  intake  tubes  which  minimize  flow  disturbance  and  enable 
concentrations  to  be  measured  to  within  1  centimetre  of  the  bed. 
This  includes  much  of  what  has  traditionally  been  regarded  as 
bed  load.  Intakes  were  situated  at  seven  different  levels  with 
a  maximum  distance  of  0.55  m  between  top  and  bottom  intakes. 

Filling  of  2.2  litre  glass  jars  via  the  intake  tubes  takes  place 
over  a  time  of  about  3.5  minutes.  The  filling  and  suction  through 
the  intakes  is  driven  simply  by  the  pressure  difference  between 
the  tops  of  the  jars  and  an  elevated  air  outlet  from  the  jars.  The 
distance  from  intakes  to  the  jars  is  1.25  m  and  the  jars  are  always 
oriented  so  as  to  minimize  flow  disturbances. 
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NEARSHORE  PROCESSES 


The  nearshore  zone,  defined  here  as  the  region  of  the  inner 
continental  shelf  over  which  shoaling  unbroken  waves  agitate 
sediment,  is  the  fundamental  'valve'  linking  shelf  and  surf  zone 
processes,  (Fig.l) .  The  nearbottom  hydrodynamic  processes  of 
this  region  determine  the  rates  and  directions  of  nearshore 
sediment  transport  while  the  configuration  of  the  nearshore 
profile  itself  determines  the  extent  of  shallow  water  wave 
modifications  and,  hence,  influences  the  amount  of  energy  which 
reaches  the  surf  zone.  We  will  confine  our  attention,  in  this 
discussion,  to  the  region  between  the  surf  zone  and  the  20  metre 
depth  contour. 

Characteristics  of  the  nearshore  profile 

Figure  11  shows  the  nearshore  profile  fronting  Eastern 
Beach.  Several  important  features  are  evident:  (1)  overall, 
the  profile  exhibits  a  pronounced  concave-upward  shape,  typical 
of  wave-built  "equilibrium"  profiles  (e.g.  Dean  1977) .  (2)  The 

innermost  200  metres  of  the  profile  consists  of  a  shallow  plateau 
2  to  3  metres  deep  which  is  surmounted  by  the  surf  zone  and  its 
associated  bar  and  trough  features.  (3)  Seaward  of  the  inshore 
plateau,  the  bed  descends  steeply  to  a  depth  of  about  15  metres 
but  with  gradients  decreasing  progressively  seaward  yielding 
an  exponential  profile  shape.  (4)  Seaward  of  a  rather  sharp 
break  in  slope  at  15  metres  (x  s  600m) ,  the  profile  flattens 
onto  a  near-horizontal  outer  nearshore  "plain";  the  20  metre 
depth  contour  lies  nearly  3  kilometres  seaward  of  the  15  metre 


contour. 
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Dean  (1977),  among  others,  has  shown  that  the  wave-dominated 
nearshore  equilibrium  profile  can  be  expressed  in  the  general  form 

h  =  a  xm  (1) 

where  h  is  local  depth,  x  is  distance  seaward  and  a  and  m  are  empi¬ 
rical  constants.  The  exponent,  m,  has  been  theoretically  evaluated 
by  Dean  (1977)  in  terms  of  two  alternative  equilibrium  conditions:  (1) 
equal  energy  dissipation  per  unit  volume  of  water,  and  (2)  equal  energy 
dissipation  per  unit  surface  area.  For  model  '1)  the  expected  value 
is  m  =  0.67  and  for  model  (2)  the  expected  value  is  m  =  0.4.  The  co¬ 

efficient,  a,  depends  on  local  environmental  conditions  and  particularly 
on  sediment  fall  velocity  w^  (grain  size)  and  dominant  wave  period,  T. 
Although  these  models  are  theoretically  strictly  applicable  to  surf 
zone  profiles.  Dean  (197’’)  has  shown  that  a  set  of  502  actual  nearshore 
profiles  from  North  America  exhibit  close  fit  to  the  model 

h  =  a  x0-67  (2) 

We  compared  the  nearshore  profile  off  Eastern  Beach  to  both  of 
Dean's  models  (1)  and  (2)  above.  The  inner  part  of  the  profile 
(h  £  15  metres)  shows  good  fit  to  either 

h  =  0.287  x067  (3) 


or 


1.233  x 


0.4 


(4) 
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Figure  11.  The  nearshore  profile  off  Eastern  Beach  based  on  surveys 
in  May,  1981.  Locations  of  current  meter  installations 
are  indicated . 
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Figure  12.  Comparison  of  the  observed  nearshore  profile  with  two 
"equilibrium"  forms  (Dean,  1977).  Note  the  abrupt 
divergence  at  h  ■  15m. 
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but  the  latter  model  (equation  4)  gives  the  best  fit.  Seaward  of 
h  =  15  metres  the  observed  profile  diverges  appreciably  from  both 
forms  as  illustrated  by  Figure  12.  This  suggests  that  the  outer 
nearshore  plain  is  morphodynamically  discontinuous  from  the  inner 
nearshore  profile.  It  appears  that  the  outer  nearshore  plain 
is  a  planation  surface  composed  of  pre-Holocene  sediments,  the  upper 
layers  of  which  are  highly  active  today  but  are  not  in  active  shore- 
normal  exchange  with  the  inner  nearshore. 

The  distribution  of  surface  sediment  characteristics  across  the 
nearshore  profile  (Table  1)  also  reveals  an  abrupt  and  dramatic  change 
in  sediment  size  just  seaward  of  the  15  metre  contour  and  supports  the 
suggestion  that  the  flat  outer  segment  of  the  nearshore  zone  is  morpho¬ 
dynamically  distinct  from  the  inner  nearshore.  From  the  10  metre  con¬ 
tour  shoreward,  the  mean  grain  size  remains  roughly  constant  at  about 
0.3-0. 4mm  (1.7-1.31?)  .  Seaward  of  15  metres  depth  grain  size  increases 
appreciably  to  over  0.65mm  at  17.5  metres  and  over  0.8  mm  at  20  metres 
(Table  1) .  The  material  seaward  of  15  metres  is  also  heavily  iron 
stained  and  numerous  whole  shells  and  cobbles  are  scattered  over  the 
bed. 

Diver  observations  at  20  metres  and  at  10  metres  revealed  large, 
symmetrical  wave-generated  ripples  at  both  positions  and  on  two 
separate  occasions.  The  ripples  had  heights  of  10  centimetres,  lengths 
(between  crests)  of  60  centimetres  and  were  very  sharp  crested.  The 
coarsest  material  (cobbles  and  whole  shells)  accumulated  in  ripple 
troughs.  The  large  ripples  and  coarse  sediment  of  the  outer  nearshore 
plain  in  combination  produced  a  very  rough  surface. 
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TABLE  1 

NEARSHORE  SAND  SIZES 


Sample 

Depth 

Mean  diameter 
(total  sample) 

Mean 

(sand 

only) 

diameter 

fraction 

Percent 

CaCO 

3 

Comments 

metres 

0  mm 

0 

mm 

0 

0 

1.13 

0.45 

1.13 

0.45 

5.0 

Beach  face 

1 

2 

1.39 

0.38 

1.43 

0.37 

3.0 

Surf  zone 

2 

3 

1.43 

0.37 

1.48 

0.36 

2.14 

Rip  bay 

3 

3 

1.67 

0.31 

1.72 

0.30 

2.24 

Outer  "Plateau" 

4 

10 

1.87 

0.27 

1.86 

0.27 

5.66 

Cobbles 

5 

10 

1.38 

0.39 

1.36 

0.39 

3.50 

present 

6 

10 

1.34 

0.40 

1.34 

0.40 

4.16 

with  some 

7 

15 

0.84 

0.56 

1.13 

0.45 

6.5 

whole  shells 

8 

17.5 

-0.16 

1.11 

0.63 

0.65 

5.56 

Cobbles  present 
with  whole  shells 

9 

20 

0.02 

0.98 

0.23 

0.85 

6.36 

II 

10 

20 

0.28 

0.82 

0.32 

0.80 

4.4 

II 

11 

20 

0.05 

0.97 

0.36 

0.78 

6.94 

It 

Observed  nearshore  currents 

The  duration  of  our  field  observations  at  20  metres  and  at  10 
metres  was  limited  by  instrument  malfunction.  The  truncation  of  re¬ 
cord  length  at  the  inner  (10  m)  station  (and  possibly  the  outer 
station  as  well)  appears  to  have  been  caused  by  the  passage  of  a 
strong  gale.  In  any  event,  we  obtained  records  of  five  15-minute 
bursts  (5  hours)  from  20  metres  and  35  bursts  (35  hours)  from  10 
metres  and  the  data  from  these  records  are  good.  From  spectral  and 
statistical  analyses  of  each  burst  we  can  draw  some  worthwhile  - 
albeit  tentative  -  inferences  about  the  magnitudes  and  directions  of 
benthic  transports  due  to  tides  and  waves.  The  locations  of  the 
h  =  20  metres  and  h  «  10  metres  current  meter  installations  are  in¬ 
dicated  in  Figure  11.  In  both  cases  the  sensors  were  located  at  1 


metre  above  the  bed 


28 


Table  2  summarizes  the  important  current  statistics  associated 
with  each  burst.  The  corresponding  deep  water  tide  level,  r)fc,  as 
observed  from  the  Kingfish  B  offshore  oil  and  gas  field  (Fig.  2)  is 
also  indicated.  The  mean  speeds  and  directions  are  the  time-averaged 
drifts  of  each  burst.  The  variance  <7^  of  speed,  V,  for  each  burst 
has  been  computed  from  the  instantaneous  (2  second  intervals)  values 
in  a  burst.  The  root  -  mean-squared  speeds  V  were  computed  from 


2  4 

v  =  <  V  > 
rms 


(5) 


(where  the  brackets  denote  time  averaging  over  a  burst) . 


The  time-averaged  means  of  the  instantaneous  speeds  cubed  <  V3>  are 
also  given. 

-+■ 

The  instantaneous  velocities,  V,  were  segregated  into  their  shore 

normal,  u  (positive  onshore)  and  shore  parallel,  v  (positive  toward 

west)  components.  The  corresponding  values  of  <u>,  u^^,  <  u3>, 

<v>,  v  ,  and  <  v3>  are  presented.  The  values  of  <V>,  <u>  and 
rms 

<v>  reflect  the  contributions  of  tidal  and  wind  driven  currents 

while  V  ,  <V3>,  u  ,  <  u3>,  v  ,  and  <  v3>  indicate  the  magnitudes 
rms  rms  rms 

of  oscillatory  flows  due  to  incident  waves  (T  =  3  -  30  secs)  and  in¬ 
fragravity  oscillations  (T  =  30  -  300  secs) . 

The  most  prominent  features  of  the  net  flows  <V>,  <u>,  <v> 
were  the  strong  dominance  of  the  shore  parallel,  <v>,  components 
(the  ratio  <|v|>/<|u|>  averaged  2.3  at  h  =  20m  and  3.5  at  h  =  10m) 
and  the  periodic  reversals  with  the  phase  of  the  semidiurnal  tide. 
These  tendencies  are  evident  from  the  stick  diagram  and  progressive 
vector  diagram  for  V  at  h  =  10m  (Fig.  13)  amd  from  the  time  series 
plots  of  <u>  and  <v>  (Fig.  14).  From  Table  2,  currents  can  be  seen 


Current  statistics  for  15  minute  bursts  from  h  =  20m  and  h  =  10m 

off  Eastern  Beach 
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Figure  13.  Vector  "stick"  diagram  and  progressive  vector  diagram 
of  burst-averaged  benthic  currents  at  h  =  10m.  Tidal 
reversals  are  apparent. 
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Figure  14 


(a)  Time  series  plots  of  burst-averaged  absolute  current 
speed,  and  the  shore-normal,  <u>,  and  shore-parallel, <v>, 
components  of  benthic  currents  at  h  *  10m;  (b)  current  rose 
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to  be  roughly  in  phase  with  tidal  elevation,  nt-  Consistent  with  the 
model  of  a  long  wave,  progressive  along  shore  from  east  to  west  (that 
is,  from  the  Tasman  Sea  into  Bass  strait;  Easton,  1970)  currents  set 
roughly  toward  the  west  at  high  tide  and  roughly  toward  the  east  at 
low  tide.  The  rotational  change  in  direction  involves  brief  periods 
of  relatively  weak  onshore-offshore  transport. 

For  the  fricti.onless  case ,  the  tidal  current  speed  V  is  given  by 
linear  long  wave  theory 

V  =  a  cos  w  t  (6) 

where  a  is  amplitude,  g  is  acceleration  of  gravity,  oj  is  radian  fre¬ 
quency  (27T/T,  T  is  period  -  in  this  case  12.42  hours)  and  t  is  time. 
For  the  tidal  amplitudes  of  0.5m  -  0.65m  which  prevailed  at  the  time 
of  the  observations,  the  frictionless  velocity  maxima  would  have  been 
approximately  0.45m  sec-1  at  h  *  20  n  and  0.64  m  sec  1  at  h  =  10m. 
However,  one  expects  rather  thick  (on  order  of  10  m)  turbulent  boun¬ 
dary  layers  for  tidal  currents  (e.g.  Komar,  1976)  and  since  the  cur¬ 
rents  were  observed  at  1  metre  above  the  bed  in  the  lower  part  of  the 
boundary  layer,  considerable  frictional  retardation  is  expected.  It 
is  thus  understandable  that  observed  currents  were  slower  than  the 
frictionless  model  predicts.  In  addition  to  the  bottom  boundary 
layer  it  also  appears  that  a  fairly  wide  lateral  nearshore  boundary 
layer  exists  to  impede  tidal  currents  with  increasing  proximity  to 
the  surf  zone.  This  is  suggested  by  the  fact  that  during  the  five 
hour  period  that  the  20  m  current  meter  was  operational,  mean  current 
speeds  were  ^oout  25  percent  weaker  at  10  m  than  at  20  m  (Table  2) . 

A  similar  shoreward  retardation  of  longshore  tidal  currents  is  des¬ 


cribed  from  a  macrotidal  environment  by  Wright  et.  al.  (in  press). 
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An  examination  of  Figures  13  and  14  and  Table  2  reveals  that  the 
easterly-setting  currents  associated  with  the  low  tide  phase  are  con¬ 
sistently  stronger  than  the  westerly  flows  which  prevail  at  high  tide. 
The  dominance  of  easterly  flows  is  particularly  obvious  from  the  cur¬ 
rent  rose  shown  in  Figure  14b  and  the  progressive  vector  diagram  shown 
in  Figure  13.  The  apparent  asymnetry  of  tidal  currents  and  the  re¬ 
sultant  net  west-to-east  transport  probably  reflects,  at  least  in  part, 
the  regional  circulation  and  wind-driven  transport  (see  discussion  in 
Fandry,  1982).  The  pronounced  increase  in  the  low-tide  easterly  flows 
and  corresponding  decrease  in  high-tide  westerly  flows  evident  in  the 
last  10  bursts  from  h  =  10  m  (1434  hrs  -  2334  hrs  on  7  May;  Table  2, 
Figs.  13  and  14)  were  related  to  southwesterly  winds  which  were  in¬ 
tensifying  over  this  period.  Since  westerly  and  southwesterly  winds 
are  dominant  for  the  region,  we  may  infer  that  the  observed  net 
currents  are  typical,  at  least  in  terms  of  direction,  of  those  which 
prevail  over  most  of  the  year. 

It  is  apparent,  from  the  V  and  maximum  speed  values  shown  in 

rms 

Table  2  that  the  instantaneous  velocities  within  a  burst,  reflecting 
the  contributions  of  wave  and  infragravity  oscillations,  dominate  over 
the  burst-averaged  speeds  and  hence  probably  play  the  dominant  roles 
in  agitating  nearshore  sediments.  Stick  diagrams  of  instantaneous 
velocities  for  a  typical  pair  of  bursts  and  time  series  plots  of  the 
associated  instantaneous  u  and  v  components  at  h  =  20  m  and  h  =  10  m 
are  shown  in  Figure  15a-c.  Histograms  indicating  the  frequency  of 
occurrence  of  different  speed  classes  within  bursts  are  shown  in 
Figure  16a  and  b  and  corresponding  power  spectra  are  shown  in  Figure 
17  a  and  b.  The  effects  of  the  waves,  causing  bidirectional  current 
oscillations,  are  readily  apparent  from  Figure  15  as  the  spectra 
(Fig.  17)  indicate,  the  dominant  energy  is  in  the  swell  wave  band 
(0.07  -  0.1  hz) .  Our  observations  as  well  as  routine  observations 


by  the  Victoria  Department  of  Public  Works  at  Lakes  Entrance  indicated 
that  the  dominant  swell  during  the  period  the  currents  were  measured 
arrived  at  roughly  45°  to  the  trend  of  the  coast.  Hence  both  u  and 
v  components  are  significant.  However,  the  interaction  of  the  oblique 
wave  oscillations  with  shore  parallel  net  tidal  currents,  apparently 
has  the  effect  of  causing  the  forward  (wave  crests)  portion  of  the 
orbital  excursion  to  be  deflected  more  onshore  while  the  reverse  flows 
(wave  troughs)  are  deflected  more  alongshore.  This  effect  is  con¬ 
spicuous  in  Figure  15a. 

To  illustrate  the  effects  of  southwesterly  gales  on  near-bottom 
oscillatory  currents,  it  is  interesting  to  compare  the  results  from 
burst  26,  recorded  shortly  before  the  passage  of  a  gale  with  the  re¬ 
sults  from  burst  30,  obtained  after  the  gale  had  become  fully  estab¬ 
lished.  Figure  18  a-c  shows  the  stick  diagrams  and  associated  u  and 
v  time  series  from  the  two  bursts  at  h  =  10m?  Figure  19  a  and  b  shows 
the  corresponding  frequency  histograms  of  speeds  and  Figure  20  a  and 
b  shows  the  comparative  power  spectra.  From  Figures  18  and  19  an  app¬ 
reciable  strengthening  of  oscillatory  flows  following  the  onset  of  the 
gale  is  obvious.  Most  importantly,  the  gale  apparently  produced  strength¬ 
ened  onshore  transport  and  easterly  transport.  The  spectra  show  that 
the  gale  caused  a  substantial  broadening  of  the  bandwidth  of  wave  os¬ 
cillations. 

From  spectral  analyses  of  the  instantaneous  u  and  v  components 
in  each  burst  (e.g.  Figs.  17  and  20  ),  it  is  possible  to  determine  the 
dominant  periods  of  the  waves  as  well  as  the  characteristic  magnitudes 
of  the  bottom  velocities.  It  is  also  possible  to  segregate  the  in¬ 
fragravity  oscillations  from  the  incident  waves  and  to  estimate  their 
corresponding  velocity  amplitudes. 
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Figure  20.  Power  spectra  of  benthic  currents,  u  and 
before  (a)  and  during  (b)  a  gale. 
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The  peak  periods,  T  ,  of  the  incident  waves  are  the  periods  at 
P 

which  the  spectral  density  (variance  per  hertz)  maxima  occur.  The 

velocity  amplitudes  ug  and  vg  corresponding  to  the  "significant"  in- 

2  2 

cident  waves  are  estimated  from  the  total  variances  a  and  0  of  u 

u  v 

and  v  in  the  frequency  band  0.033  -  0.25  hz  (T  =  4  -  30  secs)  from 

u  =  2o  and  v  =  2a  , 
s  u  s  v 

(7) 

(4s  <  T  <  30s) 

The  amplitudes  of  the  infragravity  oscillations  are  similarly  estimated 
from  the  total  variance  corresponding  to  periods  greater  than  30 
seconds.  Since  the  root  mean  squared  values  (u  ,  v  )  are  defined 


u 

rms 


/T  a 


v 

rms 


/T  a 

u 


u  and  v  are  related  to  u  and  v 
s  s  rms  rms 


by 


(8) 


u  =  /T  U  ,  V  =  /Tv 
s  rms  s  ms 


(9) 


The  total  velocity  amplitude,  Vg,  is  thus 

V  *  (u2  +  V2)*3  (10) 

s  s  s 

Knowing  the  velocity  amplitude,  V  ,  wave  period,  T  ,  and  local 

s  P 

depth,  h,  we  can  estimate  the  significant  wave  heights  Hg  from  the 
observed  currents  by  using  linear  wave  theory 
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H  =  V  T  sinh  kh/v 
s  s  p 


where  k  is  the  wave  number  i2n/L,  where  L  is  wave  length)  which  is 

related  to  the  radian  frequency  u)  (=  2tt/T  )  by  the  dispersion  rela- 

P 


k  =  0)  /g  tanh  kh 


or  for  deep  water 


=  wVg 


and  kh  may  be  approximated  from  the  explicit  expression  following  Nielsen 
(in  press) 


kh  —  /k  h  Cl  +  k  h  +  4=*-  (k  h)^J 


The  energy  density,  E,  corresponding  to  the  waves  of  height,  H  ,  is 

s 


E  =  ir  p  g 


(p  is  water  density)  and  the  energy  flux,  P,  is 


P  =  E  C 


where  is  group  velocity  which  is  related  to  phase  speed,  C  =  u>/k, 
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C 

g 


(17) 


Table  3  lists  observed  T  ,  ug,  vg  and  Vg  values  together  with 
estimated  H  and  P.  Deep  water  wave  conditions,  H  ,  T  ,  as  recorded 

g  oo  co 

from  Kingfish  B  are  also  presented.  The  components  u^n,  vj_n*  vin 
and  estimated  infragravity  height,  were  computed  from  the  variance 

in  the  infragravity  band  (30s  <  T  <  300s).  From  the  last  column 
in  Table  3  it  can  be  seen  that  infragravity  oscillations  over  the 
nearshore  zone  were  roughly  20%  to  30%  as  high  as  the  incident  waves. 
The  column  labeled  £  indicates  the  peakedness  of  the  wave  oscillatory 
currents;  it  is  given  by 


(18) 


For  a  pure  sinusoid  ?  =  1.20;  larger  values  express  the  irregularity  of 
the  waves,  due  either  to  shallow  water  deformation,  local  wind  effects 
or  both.  This  parameter  will  be  deployed  later  in  the  discussion  of 
wave  modification. 

The  information  and  discussion  just  presented  support  the 
following  general  conclusions  regarding  the  benthic  flow  field  of 
the  nearshore  zone : 

(1)  The  strongest  flows  are  the  oscillatory  currents  related  to  the 
shoaling  wind  waves  and  swell. 

(2)  Shore  parallel  reversing  tidal  currents  are  second  in  importance 
to  the  oscillatory  wave-induced  flows  in  contributing  to  the 
nearshore  current  regime.  Flows  toward  the  east  accompanying 
the  low  tide  phase  tend  to  exceed  flows  toward  the  west  which 
accompany  the  high  tide  phase. 

(3)  Over  the  observation  period,  there  was  a  net  transport  toward 


the  east. 


Wave  statistics  for  each  burst  from  h  =  20m  and  h  =  10m;  see  text  for  explanation 
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(4)  Quasi-oscillatory  flows  at  infragravity  frequencies  are  of 
tertiary  importance. 

Implications  for  sediment  transport 

In  order  for  benthic  flows  to  cause  sediment  transport,  the 
shear  stresses,  t  associated  with  these  flows  must  overcome  the 
stabilizing  force  of  gravity  on  the  sediment  particles.  An  im¬ 
portant  parameter  for  expressing  the  shear  stress  relative  to  the 
immersed  weight  of  the  particle  is  the  dimensionless  Shields  para¬ 
meter  9  (Shields,  1936?  Komar  1976)  crlven  by 

0  =  T/(p  -p)gD  (19) 

s  s 

where  p  and  p  are  respectively  sediment  and  water  densities  and  D 
s  s 

is  sediment  particle  diameter.  For  wave  induced  oscillatory  flows 
the  shear  stress  is 

x  =  h  pf  v2  =  *jpf  (a  cu) 2  (20) 

wo  w  s 

where  f  is  a  friction  factor,  V  is  orbital  velocity  maximum  and 
w  o 

ag  is  orbital  semi -excursion  length  at  the  bed 

a  =  H/(2  sinh  kh)  (21) 

s 

Historically  (and  conventionally)  the  friction  factor,  fw,  has  been 
assumed  constant  with  a  value  around  0.02  (e.g.  Bretschneider  and 
Reid,  1954).  However,  measurements  by  Carstens  et.  al .  (1969)  and 
Lofquist  (1980)  have  shown  that  f  can  have  a  large  range  of  values 
(up  to  0.3)  depending  on  the  length  scale  of  bed  roughness,  X, 
relative  to  the  orbital  semi -excursion  ag  (see  appendix  B  for  dis- 
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cussion).  Swart  (1974)  gives 

f  =  exp  [5.213  (X/a  )  '19“  -  5.977]  (22) 

w  s 

For  a  fixed  (immobile) ,  unrippled  bed  the  roughness  length  is  about 

X  =  2.5  D  (23) 

s 

(Nielsen,  in  press)  which,  for  D  ^  1  ran  (X  =  .0025  m)  and  a  =  0.20m 
(H  =  .5m,  T  =  10  sec,  h  =  20  m)  the  friction  factor,  estimated  from 
equation  (22)  would  have  a  "conventional"  value  on  the  order  of  0.02. 
For  initiation  of  sediment  motion  by  waves,  it  is  normal  to  use  the 
skin  friction  Shields  parameter,  0‘,  obtained  by  substituting  (20) 
into  (19)  and  using  (23)  for  roughness  (e.g.  Nielsen  1979). 

Numerous  studies  of  the  critical  or  threshold  value,  0  ,  which 

c 

must  be  exceeded  before  sediment  transport  can  be  initiated  suggest 
that  the  critical  values  are  about  0  =  0.04  -  0.06  for  sand  trans¬ 
port  by  both  steady  and  oscillatory  flows  over  smooth  beds  (Sternberg, 
1972;  Komar,  1976;  Komar  and  Miller,  1975)  over  rippled  beds,  the 
critical  value  may  be  as  low  as  8^  ~  0.03. 

For  steady  or  quasi-steady  flows  (e.g.  tidal  currents),  Sternberg 
(1972)  found  that  the  drag  coefficient,  relating  current  speed, 

<V>^qq  at  1  metre  above  the  bed  to  shear  stress,  X,  has  a  mean  value 
of  3.1  x  10“ 3  or  in  other  wcrds  for  steady  (non-wave)  currents 

X  a  3.1  x  10"3  p<V>2  (24) 

Substituting  observed  burst- averaged  current  speeds,  <V>,  (Table  2) 

into  the  above  relationship  for  X,  using  observed  sediment  sizes 

(Table  1)  and  assuming  p  =  2650  kg  m-3  we  estimate  the  associated 

s 
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for  D  *  0.2  mm  and  9  »  0.008  for  D  =  0.8  mm.  Both  values  are 

s  s 

subcritical  and  it  may  be  inferred  that,  if  we  accept  Sternberg's 
(1972)  estimate  of  the  drag  coefficient,  the  observed  non-wave-in- 
duced  flows  are  unable  to  initiate  significant  nearshore  sand  trans¬ 
port. 

As  concluded  previously,  the  oscillatory  velocities,  V  ,  due  to 
the  waves  are  dominant  over  the  nearshore  profile  and  we  may  expect 
these  motions  to  play  the  dominant  role  in  agitating  nearshore  sedi¬ 
ments.  We  will  evaluate  the  likelihood  of  sediment  being  entrained 
by  waves  using  the  semi-empirical  relationships  proposed  by  Komar  and 
Miller  (1975) .  Komar  and  Miller  (1975)  used  a  Shields-type  criterion 
or  mobility  number  of  the  form 


(pg  -  p)  g  ds 


(26) 


where  Vq  is  maximum  orbital  velocity  (approximated  in  our  case  by  V  ) . 
They  found  that  for  turbulent  boundary  layers  on  an  un rippled  bed 
the  critical  value,  ¥  ,  necessary  to  initiate  motion  is 


V 

c 


0.46 IT  (d./D  ) 
0  s 


k 


(27) 


where  dQ  is  the  bottom  water  particle  excursion  length,  dQ  =  2ag 
(eg.  21) . 

On  the  basis  of  the  above  criterion,  ¥  ,  bottom  orbital 

c 

velocities  of  0.18  m  sec-1  would  be  sufficient  to  move  the  finer 
fractions  (Dg  =  0.2  mm).  This  velocity  was  exceeded  by  nearly  half 


h 
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of  the  observed  V  values  shown  in  Table  3.  We  can  estimate  that 
s 

10  second  waves  would  require  a  local  height  of  0.7  m  at  h  =  20  m 
and  0.42  m  at  h  =  10  m  in  order  for  the  associated  orbital  velocities 
to  exceed  0.18  m  sec-1.  For  the  much  coarser  mean  grain  sizes  of 

0.8  nan  at  h  =  20  m  and  0.4  mm  at  h  =  10  m  the  respective  "threshold 

velocities"  would  be  0.33  m  sec-1  and  0.25  m  sec-1;  the  corresponding 
"threshold"  wave  heights  (T  =  10  sec)  are  H  *'  1.3  u  at  h  *  20  n  and 

H  =  0.58  m  at  h  =  10  m.  The  4*  criterion  as  expressed  by  equation 

27  assumes  a  smooth  bed.  The  nearshore  bed  off  Eastern  Beach  was 
observed  to  be  extremely  rough  with  sharp  crested  wave  ripples  10 
centimetres  high.  Over  a  rippled  bed  sediment  entrainment  should  take 
place  at  lower  orbital  velocities.  Diver  observations  confirmed  that 
oscillatory  motions  were  entraining  sediment  clouds  over  ripples. 

We  conclude  that  waves  are  probably  able  to  agitate  nearshore 
sediments  shoreward  of  the  20  metre  contour  for  most  of  the  time. 
However,  agitation  of  sediments  by  oscillatory  currents  does  not 
necessarily  imply  any  net  transport.  We  are  concerned  now  with 
inferring  the  likely  direction  of  net  nearshore  transport.  We  will 
not  attempt  to  estimate  the  absolute  rates  of  transport  because  the 
data  are  insufficient  to  allow  this. 

Many  recent  studies  of  nearshore  sediment  transport  attempt  to 
model  the  transport  rate,  q,  in  terms  of  modified  forms  of  Bagnold's 
(1963)  energetics  approach  (e.g.  Sternberg,  1972;  Seymour,  1980; 
Bailard,  1981;  Bailard.  and  Inman,  1981).  This  approach  considers  q 
to  be  proportional  to  the  bed  "fluid  power"  that  is 

q  oc  x  V  (28) 

This  approach  has  been  applied  most  often  to  bedload  transport,  qfc, 
although  Bailard  (1981)  has  proposed  a  "total  load"  version.  If  the 
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net  tidal  and  wind  driven  currents  operated  alone  without  wave  effects, 
then  these  currents  would  be  the  source  of  the  shear  stresses  and  we 
would  expect 


qfa  <*  <  V3>  (29) 

However,  we  have  shown  that  waves  are  the  dominant  source  of 
bed  shear  stress.  Prediction  of  net  sediment  transport  by  a  com¬ 
bination  of  waves  and  currents  is  a  very  complicated  matter  and  at 
present  there  is  no  model  in  existence  which  at  the  same  time  is 
both  theoretically  based  and  in  agreement  with  experimental  evidence. 

As  early  as  1963,  Inman  and  Bowen  (1963)  showed  that  the  net 
sediment  transport  under  a  combination  of  waves  and  currents  will  just 
as  often  be  against  the  net  current  as  with  it.  This  means  that  no 
transport  model  by  which  the  net  transport  is  proportional  to  the  net 
flow  will  be  in  general  agreement  with  the  empirical  evidence. 

More  such  evidence  has  since  been  obtained  by  many  detailed 
flume  experiments  at  the  Technical  University  in  Delft;  see  van  de 
Graaff  (1980) .  This  effect  is  due  to  the  fact  that  the  transport  con¬ 
tribution  from  oscillatory  velocity  and  concentration  components 

U  and  C  will  often  dominate  over  the  contribution  of  the  steady  com- 
n  n  1 

ponents  , 

Q  =  C(z)<u>(z)  dz  (30) 

o  o 

This  dominance  is  due  to  the  fact  that  the  oscillatory  velocity  com¬ 
ponents  will  have  a  much  thinner  boundary  layer  than  the  steady  flow, 
and  will  thus  penetrate  to  the  lowest  elevations,  where  the  concen¬ 
trations  are  largest. 


k. 


> 


i 

i 


The  direction  of 


t+T  h 

i  dt  |  Cn  (z,t)  (z,t)  dz  (31) 

t  o 

will  depend  on  the  phase  relations  between  the  simple  harmonics 
(z,t)  and  Un(z,t)  and  as  yet  there  is  no  model  which  is  capable 
of  predicting  these  phase  relations  on  theoretical  grounds. 

The  model  presented  by  Nielsen  et.  al.  (1978)  and  Nielsen  (1979) 
predicts  these  phase  relations  and  predicts  net  sediment  fluxes  that 
are  in  agreement  with  the  Dutch  measurements,  with  appropriate  choice 
of  parameters.  However,  this  model  is  based  on  the  diffusion  equation 
and  recent  empirical  and  theoretical  work  shows  that  the  entrainment 
of  sand  by  a  wave  dominated  flow  is  net  a  diffusion  process  (Nielsen 
in  prep. ) . 

In  the  direction  perpendicular  to  the  wave  motion  the  problem 
is  not  so  conplex  because  there  are  no  oscillatory  velocity  components 
and  therefore  the  transport  is  given  reasonably  by 

h 

c  (z)  <V>  (z)  dz  (32) 

o 

Even  here  however  c  will  be  determined  by  the  combined  shear  stress 
from  the  waves  and  currents  and  there  is  at  present  no  theory  which 
gives  a  satisfactory  description  of  the  boundary  layer  structure  for 
combinations  of  waves  and  currents.  The  existing  models  e.g.  Madsen 
and  Grant  1979  all  assume  a  constant  or  linearly  increasing  eddy 
viscosity  which  is  in  clear  disagreement  with  the  eddy  viscosity  mea¬ 
surements  referred  to  by  Lundgren  (1972)  and  with  the  velocity  fluctuation 
measurements  by  Nakato  et.  al.  (1977)  and  Macdonald  (1977) .  The 
latter  two  both  found  that  the  intensity  of  vertical  velocity  fluc¬ 


tuations  decreased  exponentially  with  the  distance  from  a  rippled  bed, 


50 


and  such  exponentially  decreasing  velocity  fluctuations  are  not  very 
likely  to  correspond  to  a  constant  or  linearly  increasing  eddy 
viscosity. 

The  key  to  a  consistent  theoretical  description  of  wave  and 
wave/c  .’rrent  boundary  layers  over  rippled  beds  seems  to  be  the  con¬ 
sideration  of  the  fact  that  both  suspended  sand  and  turbulence 
originate  at  the  bed  and  move  upward  by  a  rhythmic  convective  process. 
The  nature  of  this  process  is  indicated  by  the  measurements  of  Nakato 
et.  al.  (1977)  and  Macdonald  (1977)  who  found  that  both  suspended 
sediment  and  turbulence  are  exponentially  distributed  with  identical 
vertical  length  scales. 

In  the  absence  of  good  theories  for  the  boundary  layer  struc¬ 
ture  of  the  combined  flow  we  are  forced  to  use  the  very  simple 
assumption  that  the  concentrations  c(z)  are  determined  by  the  waves 
alone.  This  assumption  seems  to  be  reasonably  well  justified  under 
most  conditions  even  in  the  presence  of  quite  strong  currents  (e.g. 
longshore  currents  as  discussed  in  the  next  section) . 

We,  therefore,  can  draw  no  conclusions  as  to  shore-normal 
transport.  We  can,  however,  make  some  inferences  at  least  concerning 
the  direction  of  shore -parallel  transport  if  we  regard  the  process 
as  involving  entrainment  of  sediment  by  wave  induced  oscillatory  cur¬ 
rents,  Vs»  and  net  advection  of  the  entrained  sediment  taking  place 
in  response  to  the  net  time  averaged  longshore  current  <v>.  An 
energetics  form  for  the  superimposition  of  wave  shear  stress  and 
steady  currents  is  given  by  Seymour  (1980)  for  the  case  of  inter¬ 
acting  waves  and  tidal  currents.  Seymour  (1980)  suggests  that 

q  a  «  v3»  +  <<v><V2  >>  (33) 

b  s 

(The  second  set  of  brackets  denotes  time  averaging  over  all  bursts) . 
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Seaward  of  the  surf  zone,  the  effect  of  wave  oscillations  over 
a  rippled  bed  is  to  suspend  sediment  in  vortices  above  the  bed  and 
there  is  considerable  field  and  laboratory  evidence  that  sediment 
transport  by  unbroken  waves  over  a  rippled  bed  takes  place  largely  in 
intermittent  suspension  {Nielsen,  1979;  in  press) .  If  the  concen¬ 
tration  of  suspended  sediment  near  the  bed  is  Cg  then  the  alongshore 
suspended  load  transport  rate  will  be 

Q  a  <  v>  C  (34) 

s  s 

Field  observations  of  suspended  sediment  concentration  profiles  mea¬ 
sured  under  waves  (Nielsen  and  Green,  in  press;  Nielsen,  1979,  see 
discussion  of  surf  zone  transport  later  in  this  report)  indicate  that 
Cs  is  roughly  proportional  to  the  cube  of  the  Shields  parameter  and, 
hence,  to  the  sixth  power  of  the  wave  orbital  velocity,  that  is 

C  «  0’3  T3  “  V6  (35) 

s  s 

We  may  therefore  consider  that 

q  «  <v  >  V6 
s  s 

Table  4  lists  the  observed  values  of  all  of  the  sediment  trans¬ 
port  indices  just  described,  time  averaged  over  the  35  hour  observation 
period  at  h  =  10  m.  All  shore  parallel  indices  are  negative  in¬ 
dicating  a  net  transport  toward  the  east. 


L 
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Table  4:  Indices  of  relative  shore-parallel  sediment  transport 
rate  a:  h  =  10  m. 


<<  v  >> 


m3sec"3 


«  V  V2>>  «  V3>>  +  «  V  V2>>  <<  V  v6>> 

s  s  s 

3-3  3  — 3  7—7 

m  sec  m  sec  m  sec 


Shore-Parallel 

(V)  - .6x10“ 3  -.73x10" 3  -. 133x10" 2  -.23x10 


v  positive  toward  west  (alongshore) 


During  the  experiment  period,  the  easterly  setting  flows  dominated 
over  the  westerly  setting  flows.  Naturally,  we  cannot  assume,  on  the 
basis  of  a  short  sample,  that  this  is  always  the  case.  We  can,  however, 
be  reasonably  confident  of  one  important  conclusion:  tidal  currents  set 
to  the  east  at  low  tide.  As  emphasized  by  Seymour  (1980)  and  subseq¬ 
uently  by  Wright  (in  press) ,  low  tide  flows  are  likely  to  transport 
more  sediment  than  high  tide  flows  of  equal  magnitude  because  wave- 
induced  shear  stresses  at  the  bed  will  be  greater  at  low  tide  when 
depth  and,  hence  sinh  kh,  are  reduced.  The  evidence  at  hand  suggests 
that  net  sediment  transport  over  the  nearshore  profile  is  toward  the 
east. 

Nearshore  modifications  of  waves 

From  the  point  of  view  of  the  surf  zone  processes  to  be  con¬ 
sidered  in  the  following  section,  the  modifications  experienced  by 
the  incident  waves  as  they  cross  the  nearshore  zone  between  deep 
water  and  the  break  point  are  of  considerable  importance.  The  most 
important  modifications  include  (1)  shoaling;  (2)  refraction;  and 
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(3)  dissipation  by  bed  fric«.xun.  Local  wind  shear  can  also  produce 

appreciable  modification.  Considering  wave  height  variations  due  to 
those  processes,  wave  height,  H,  at  any  given  position  and  depth  re¬ 
lative  to  deepwater  wave  height,  H  is 


=  K  K  K, 
H  S  r  f 


where  are  respectively  the  shoaling,  refraction  and  friction 

coefficients. 


The  shoaling  coefficient,  K  is 

s 


K  -JtP-  = 
S  VC 

g 


0.5  C 


where  n  =  ^[l  +  sinh~2k'ti  2111(3 '  the  case  of  straight,  parallel  depth 
contours  (a  condition  approximately  met  at  the  field  site)  the  re¬ 
fraction  coefficient  is 


where  a  and  are  the  local  and  deepwater  angles  of  wave  incidence 
which  vary  in  accordance  with  Snell's  law 


=  constant 


Using  the  explicit  formulae  developed  by  Nielsen  (in  press)  we  can 
approximate  Kg  and  from 


K  =  (8irh/L  )  exp  (trh/2L  ) 


and 
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K  =  /cosa  [l-27Th/L  (1-h/L  )2  sin2a  ] 

j r  OO  00  00  OO 


where  the  deepwater  wave  length,  L  is 


gT  /2t r 


If  we  could  imagine  the  incident  waves  to  be  unaffected  by  friction 

(which  we  cannot)  then  we  could  expect  the  "modal"  deepwater  waves 

ofHoo=3m,  T  =  8.5  sec,  =  45°  (SW)  to  break  at  an  angle  of 

16.8°  with  a  breaker  height  of  slightly  less  than  3  m  corresponding 

to  K  =  1.15  and  K  =  0.86. 
s  r 

* 

As  suggested  earlier,  the  effects  of  frictional  dissipation  over 
the  wide,  rough  nearshore  bed  are  probably  large.  Whereas  the  effects 
of  refraction  and  shoaling  alone  involve  no  change  in  energy  flux,  P, 
frictional  dissipation  results  in  a  loss  of  energy  flux  at  a  rate 


1Z  =  _  x  V  z  ->ipf  jv3|=->5pf  ^(V  ) 3  - - —  pf  £V3  (43) 

d  s  w  w  rms  w  ° 


where  s  is  horizontal  distance  along  the  direction  of  wave  propoga- 
tion  ,  V  is  the  instantaneous  bottom  orbital  velocity  and  Vq  is 
the  maximum  bottom  orbital  velocity 


o  T  sinh  kh 


In  our  preceding  discussion  of  sediment  transport,  we  estimated 
f^  from  equation  (22)  using  a  roughness  length  of  X  =  2.5  Dg  for  pur¬ 
poses  of  calculating  the  skin  friction  Shields  parameter  Q't  appro¬ 
priate  to  the  grains.  However,  in  order  to  estimate  frictional  dissi¬ 


pation,  we  must  consider  the  total  combined  roughness  for  mobile  (as 
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opposed  to  fixed)  beds  and  the  presence  of  the  ripples.  The  roughness 
X,  and  consequently  f^,  are  considerably  increased  over  mobile  rippled 
beds.  On  the  basis  of  several  sets  of  available  data  (e.g.  Carstens 
et.  al.  1968;  Lofquist  1980),  we  find  the  best  empirical  fit  for  X 
as  a  function  of  ripple  dimensions  and  grain  size  to  be 

X  =  Br\2/L  +  190/9''  -.05'  D  (45) 

r  r  s 

where  is  ripple  height,  is  ripple  chord  length  (spacing  between 
ripple  crests) ,  Q'  is  the  skin  friction  Shields  parameter.  The 
rationale  behind  this  expression  is  explained  in  Appendix  B.  Using 
observed  ripple  dimensions  off  Eastern  Beach  r)  =  0.10  m,  =  0.6  m, 

X  is  found  to  be  about  0.13  m  under  low  wave  energy  (  H  »  0.5  m)  when 
9',  and  hence  the  moving  sand  contributions  are  small.  Under  larger 
waves,  when  9'  at  h  *  20  m  is  on  the  order  of  0.20  for  the  coarser 
fractions  and  up  to  0.80  for  the  finer  fractions,  the  moving  sand 
contribution  would  increase  X  to  the  order  of  0.20m.  Substituting 
our  new  roughness  estimates  into  Swart’s  (1974)  formula  (eg.  22) 
we  obtain  estimates  of  f  ranging  from  f^  =  0.18  to  0.31;  a  typical 
rippled  bed  value  is  f  =  0.25  which  we  will  assume  to  be  roughly 
constant  over  the  nearshore  profile. 

Since  f  probably  has  a  large  value  of  ~  0.25  and  since  the 

modal  southwesterly  waves  have  a  large  initial  angle  of  incidence 

(a  't  45°)  and  hence  must  travel  large  distances  over  the  low  qra- 
00 

dient  nearshore  profile,  considerable  loss  of  energy  must  take  place 
before  the  waves  reach  the  break  point.  To  estimate  the  likely 
breaker  heights  corresponding  to  the  modal  deepwater  waves  we  must 
consider  the  combined  effects  of  refraction,  shoaling,  and  fric¬ 
tion.  We  can  do  this  bv  comouti no  chancres  over  consecutive  segments 
of  the  profile.  If  h^  and  h2  are  depths  at  the  seaward  and  landward 
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ends  of  a  segment  and  and  H are  the  corresponding  heights 
then  we  cam  approximate  the  ratio  from 


where  subscripts  1  and  2  denote  the  values  at  deeper  and  shallower 
ends  of  the  segment,  the  subscript  “  denotes  the  deepwater  value, 
dh/dx  is  bed  slope,  £  is  the  peakedness  function  (eq.  18)  which  in 
our  case  is  about  2.05  (on  average)  and  the  integral  1  is 

ci.5 

- - —  „  -  . .  (cos  a)”2-5  d  k  h  (47) 

CLS  sinh3  kh 

g 


This  integral  can  be  approximated  by 


where 


0.7  tan2^ 

(l+tan2a.)k  h, (1  -  i  k  h.) 
1  00  1  3  00  1 


(49) 


and 
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4  tan4a 

U  =  - - -  _  T  5  (50) 

(l+tan2^)2^^)2 

The  cumulative  error  of  the  integral  I  can  be  inferred  from 
the  error  on  the  corresponding  explicit  approximation  to  the  inte¬ 
grand,  and  a  detailed  error  analysis  can  be  found  in  Nielsen  (in 
press).  For  <  45°  the  induced  error  on  the  predicted  wave 
height  will  be  insignificant. 

Solving  equation  46  for  the  estimated  modal  wave  conditions  across 

successive  increments  of  the  nearshore  profile  having  constant 

gradients,  using  f  =0.25  and  the  observed  C  value  of  2.05  and 
w 

assuming  negligible  modifications  seaward  of  h  =  20  m  we  estimate 
the  modal  breaker  height  to  be  about  1  metre  and  the  breaker  angle 
would  be  about  9.6°.  Figure  21  shows  the  predicted  change  in 
heights  across  the  profile  for  shoaling  and  refraction  alone  and  for 
shoaling,  refraction  and  friction.  It  is  apparent  from  Figure  21 
that  most  of  the  wave  height  reduction  takes  place  over  the  wide,  low 
gradient  segment  netween  h  =  20  m  and  h  =  15  m. 

The  predicted  reduction  in  wave  height  to  about  one  third  the 
deepwater  value  before  the  waves  break  may  seem,  to  some,  to  be  ex¬ 
cessive.  Perhaps  we  have  overestimated  f  ?  He  can  check  this  to 

w 

a  limited  degree  because,  for  a  short  period  of  five  hours  we  have 

roughly  simultaneous  measurements  of  V  and  hence  estimates  of  H 

s  s 

from  h  =  20  m  and  h  =  10  m  (Table  3).  From  the  H  values  in  Table 

s 

3  we  can  estimate  the  dissipation  factor,  Fg,  (which  includes  f^ 
plus  any  other  dissipating  effects  which  may  be  present)  from 
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h  (m) 


Variation  of  local  wave  height,  H,  with  depth,  h, 
as  predicted  for  a  deepwater  wave  with  «  3.0, 

T  *  8.5  sec  on  the  basis  of  shoaling  and  refraction 
alone  and  shoaling  refraction  and  frictional  dissi¬ 
pation  over  the  observed  nearshore  profile. 


Figure  21 
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dh 

ds 


kJ*ll 


1 


(51) 


where  and  must  be  the  rms  values  which  are  related  to  by 

H  =  H  //T!  For  irregular  waves,  this  eauation  would  result  in 
rms  s 


overestimating  Fg  unless  we  include  the  observed  values  of  £  (in- 

a/T 

stead  of  the  "simple  harmonic  value" 


(52) 


For  the  observed  wave  conditions,  the  angle  of  incidence 
was  at  both  stations  close  to  45°  so  that  no  refraction  had  taken 
place(since  ug  and  vs  were  nearly  equal)  so  dh/ds  =  dh/dx  For 

the  five  hourly  bursts  during  which  both  inner  and  outer  current 

meters  were  recording  simultaneously  the  estimates  of  Fe  are  given  in 
Table  5 . 


»' 


Table  5. 

Burst 
F 

e 

The  dissipation  factors  obtained  in  this  way  all  substantially 

exceed  the  predicted  f  value  of  0,25.  Some  of  this  apparent  excess 

w 

dissipation  is  probably  not  real  but  is  only  a  manifestation  of  the 
non-linear  effects  of  the  shoaling  waves  (R.T.  Guza  Pers.  Com.; 
Lighthill  1978) .  Some  may  be  due  to  extra  dissipation  caused  by 
interaction  between  the  waves  and  the  tidal  currents;  and  some  may 


F  values  estimated  from  observed  current  statistics, 
e 


1 

.71 


2 

.65 


.31  .59  .52  F  =  0.56 

e 


be  due  to  the  fact  that  we  underestimated  rather  than  overestimated 


! 
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the  friction  factor,  f  . 

w 

The  effects  of  frictional  dissipation  are  non-linear;  the  rate 
of  dissipation  is  proportion^  to  the  cube  of  the  wave  height  (eq.  43). 
Hence  the  large  waves  will  experience  more  attenuation  than  the  small 
waves.  An  important  question  of  fundamental  relevance  to  beach  and 
surf  zone  variability  thus  arises:  how  does  the  nearshore  modification 
of  waves  by  shoaling,  refraction,  and  friction  affect  the  annual  dis¬ 
tribution  (temporal  variability)  of  nearshore  (or  breaker)  heights? 

The  effects  of  shoaling  and  refraction  are  linear  and  thus  only  pro¬ 
duce  shifting  and  stretching  of  the  distribution.  The  addition  of 
frictional  effects,  however  can  be  expected  to  cause  narrowing  and 
skewing  of  the  wave  height  distribution  near  the  break  (or  in  shal¬ 
low  water)  and  reduce  the  frequency  of  occurrence  of  larger  heights. 
This  effect  is  very  pronounced  when  f^  is  large  and  when  the  near¬ 
shore  profile  is  flat. 

Equation  46  gives  the  wave  height,  H2,  at  the  shallow  limit 
of  the  profile  segment  relative  to  the  height,  ,  at  the  deep  . 
(offshore)  limit  of  the  segment  after  modification  by  refraction, 
shoaling  and  friction.  For  simplicity,  we  will  rewrite  equation 
46  as 


h2  ■  G<K1>  ■  1  !  n1  h  l53> 


and  I  is  given  by  equation  47. 
Now,  if  the  distribution  of  H  is 


P  (H1  _<  H}  =  F^  (H) 


(56) 


(where  P  (h^  _<  H}  is  the  annual  probability  chat  is  less  than  or 
equal  to  a  given  height) 


and  H2  =  G(HX) 


then  the  distribution  of  H2  is 

P  (h2  <  h}  =  F2(H)  =  F1  Eg"1  (H)] 

From  (53) 


H1  "  &  ~  n  h2 


(57) 


(58) 


(59) 


SO 


-  F^Th  (60) 

and  the  distribution  of  H2  is 

P  <H2  <H>  -  P2  (H)  -  FjIjHVh’  <61) 

By  differentiating  (61)  with  respect  to  H  we  find  the  probability 
2  ' 


density  of  H 
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P2  (H) 


(ip  -  OH)* 


(62) 


where  p^  and  p2  are  the  probability  densities  of  heights  at  h^  and  h2> 

The  distribution  of  deepwater  wave  heights  in  Gippsland  is  broad.  If 
(defined  for  h.^  ®  20  metres)  is  assumed  to  follow  a  Rayleigh  distri¬ 
bution  (which  it  probably  does  not  really  do),  then. 


Fx  (H)  *  1  -  exp  [-(H/H)2] 


(63) 


where  H  is  the  mean  height  (H  =  2  m  in  our  case) . 

Then,  from  (61) 


H2 

f2  (h)  «  1  -  exp  c-  Cip  _-a-Hf2  a2  ] 


(64) 


and  from  (62)  the  probability  density  is 
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2  ip  H 

Op  -  ft  H)3  H2 


exp 


H2 

Op  -  ft  H)  2  H2 


(65) 


To  obtain  approximate  extimates  of  p2 (H)  we  use  observed  para¬ 
meter  values  for  the  Eastern  Beach  case  of  h^  *  20  m,  h2  *  2  m,  H  *  2  m, 
ip  »  1.4,  ft  «  0.59  corresponding  to  the  predicted  f^  of.  0.25  and  ft 
*  1.32  corresponding  to  the  observed  Ffi  value  of  0.56. 

Figure  22  shows  the  estimated  probability  density  p^(H)for  h^  »*  20  m 

("deepwater")  and  p2  (H)  estimated  for  the  "no  friction"  case 

(f  ,ftl  -  ( zero, zero) ,  the  "estimated  dissipation  case"  (f  ,ft>  ■  (0.25,  0.59) 
w  w 


Figure  22.  Probability  density  of  deepwater  wave  heights,  Pj (H) , 

(assuming  a  Rayleigh  distribution)  and  breaker  heights, 
PjOJ)  after  frictional  dissipation  with  f  *  0.25  or 
with  F@  ■  0.56.  Breaker  height  distribution  on  the 
basis  of  no  dissipation  is  also  indicated. 
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and  the  "measured  dissipation  case"  (fw,Q)  =  (0.56,  1.32).  The  deep¬ 
water  distribution  of  is  very  broad,  the  temporal  variability  is 
large  and  large  wave  heights  are  relatively  frequent.  If  frictional 
dissipation  did  not  operate,  the  distribution  of  breaker  heights 
would  be  shifted  somewhat  but  would  remain  broad  with  large  maxima. 

However,  a  bottom  friction  factor,  f  of  0.25  combined  with  the  low 

w 

gradient  nearshore  profile  substantially  narrows  the  distribution  of 
breaker  heights  and  reduces  the  occurrence  of  large  breakers  despite 
the  fact  that  much  larger  waves  may  occur  offshore.  For  Fg  values  of 
0.56  the  distribution  would  be  narrower  still  and  the  maximum  height 
would  rarely  exceed  1  metre. 

Whatever  the  absolute  p2 (H)  values  may  really  be,  one  qualitative 
point  is  very  clear:  the  temporal  variability  of  breaker  height  is  much 
lower  than  the  temporal  variability  of  deepwater  wave  height  and  the 
maximum  breaker  height  is  apparently  limited  to  about  2  metres  for  most 
of  the  time.  This  is  a  consequence  of  the  rough,  low-gradient  near¬ 
shore  profile  fronting  the  Gippsland  coast.  Over  steeper,  narrower  near¬ 
shore  profiles  such  as  that  fronting  the  coast  of  New  South  Wales  (e.g. 
Wright  1976)  frictional  attenuation  is  less  and  the  variability  of 
breaker  heights  is  much  broader.  This  is  reflected  in  the  stronger 
seasonal  variability  for  most  New  South  Wales  beaches  compared  to 
Eastern  Beach. 
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SURF  ZONE  AND  BEACH  PROCESSES 

The  surf  zone  and  beach  at  Eastern  Beach  occupy  only  the  innermost 
200  metres  of  the  Coastal  Boundary  Layer  (Fig.l).  However,  it  is  within 
this  comparatively  narrow  region  that  the  previously  modified  waves  ulti¬ 
mately  exoerience  their  most  rapid  dissipation,  perform  work  at  the  most 
rapid  rate,  and  produce  the  most  dramatic  and  rapid  morphologic  changes. 

After  the  incident  waves  break  at  the  outer  limits  of  the  surf  zone,  their 
energy  is  transferred  to  other  forms  of  motion  including  lower  frequency 
standing  waves  and  edge  waves,  longshore  currents,  and  rip  circulations. 

Morphodunamia  state  of  the  surf  zone  and  beach 

Figure  23  shows  the  six  most  common  morphodynamic  states  recognized 
by  Wright  and  Short  (in  press) .  The  general  model  shown  in  Figure  23  is  a 
modified  and  updated  version  of  the  models  originally  discussed  by  Wright 
et  al  (1979  a  and  b)  and  Short  (1979  a  and  b) .  The  two  extremes  shown  in 
Ficrure  23  a  and  f  are  respectively  fully  dissipative  and  highly  reflective. 
Dissipative  versus  reflective  conditions  are  distinguished  on  the  basis  of 
a  surf  scaling  parameter  proposed  by  Guza  and  Inman  (1975)  of  the  form 

e  =  a  w2/ g  tan2  5  (66) 

where  a  was  runup  amplitude  (in  Guza  and  Inman's  1975  original  equation), 
w  is  incident  wave  radian  frequency  (w  =  27T/T,  T  is  period)  g  is  accelera¬ 
tion  of  gravity  and  8  is  beach  or  surf  zone  gradient.  Wright  et  al  (1979) 
replaced  runup  amplitude  a  with  breaker  amplitude,  a^  to  differentiate  between 
reflective  and  dissipative  surf  zones  and  we  follow  that  approach  here. 

Guza  and  Inman  (1975)  and  Guza  and  Bowen  (1977)  showed  that  strong  reflection 
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Figure  23.  Morphodynamic  states  of  surf  2ones  and  beaches 

(a)  dissipative,  (b-e)  intermediate;  (f)  reflective 


SURF  ZONE  AND  BEACH  PROCESSES 


The  surf  zone  and  beach  at  Eastern  Beach  occupy  only  the  innermost 
200  metres  of  the  Coastal  Boundary  Layer  (Fiq.l).  However,  it  is  within 
this  comparatively  narrow  region  that  the  previously  modified  waves  ulti¬ 
mately  experience  their  most  rapid  dissipation,  perform  work  at  the  most 
rapid  rate,  and  produce  the  most  dramatic  and  rapid  morphologic  changes. 

After  the  incident  waves  break  at  the  outer  limits  of  the  surf  zone,  their 
energy  is  transferred  to  other  forms  of  motion  including  lower  frequency 
standing  waves  and  edge  waves,  longshore  currents,  and  rip  circulations. 

MorphodynanrLa  state  of  the  sitrf  zone  and  beaeh 

Figure  23  shows  the  six  most  common  morphodynamic  states  recognized 
by  Wright  and  Short  (in  press) .  The  general  model  shown  in  Figure  23  is  a 
modified  and  updated  version  of  the  models  originally  discussed  by  Wright 
et  al  (1979  a  and  b)  and  Short  (1979  a  and  b) .  The  two  extremes  shown  in 
Fiaure  23  a  and  f  are  respectively  fully  dissipative  and  highly  reflective. 
Dissipative  versus  reflective  conditions  are  distinguished  on  the  basis  of 
a  surf  scaling  parameter  proposed  by  Guza  and  Inman  (1975)  of  the  form 

e  =  a  w2/g  tan2  8  (66) 

where  a  was  runup  amplitude  (in  Guza  and  Inman's  1975  original  equation)  , 
a)  is  incident  wave  radian  frequency  [u  =  2g/T,  T  is  period)  g  is  accelera¬ 
tion  of  gravity  and  8  is  beach  or  surf  zone  gradient.  Wright  et  al  (1979) 
replaced  runup  amplitude  a  with  breaker  amplitude,  a^  to  differentiate  between 
reflective  and  dissipative  surf  zones  and  we  follow  that  approach  here. 

Guza  and  Inman  (1975)  and  Guza  and  Bowen  (1977)  showed  that  strong  reflection 


67 


and  associated  standing  wave  motion,  particularly  at  subharmonic  frequencies, 
and  surging  breakers  with  high  runup  prevail  when  £  <  2. 0-2. 5.  Wright 
at  al  (1979  a  and  b)  and  Wright  (1982)  found  that  highly  reflective 
natural  beaches  in  Australia  consistently  exhibit  £  values  less  than 
2. 5-3.0;  1. 5-2.0  is  typical.  When  e  >  2.5,  waves  begin  to  plunge  and 

when  e  >  20  spilling  breakers  occur  (Galvin,  1972) .  Surf  zone  width  and 
turbulent  dissipation  of  incident  wave  energy  increase  with  increasing  e . 

In  the  fully  dissipative  extreme,  high  e  values  (e  =  30-500)  prevail  across 
the  entire  surf  zone.  Breakers  are  exclusively  spilling  and  decay  progress¬ 
ively  across  the  flat  bed  to  become  very  small  by  the  time  they  reach  the 
subaerial  beach  (e.g.  Wright  et  al  1982). 

Each  of  the  intermediate  states  shown  in  Figure  23  b-e  possesses 
both  dissipative  and  reflective  elements  and  e  varies  appreciably  either 
across  the  profile  or  alongshore  (e.g.  Wright  and  Short,  in  press) .  Our 
observations  of  Eastern  Beach  indicate  that  this  beach  typically  alternates 
between  the  longshore  bar-trough  state  (Fig. 23b)  and  the  rhythmic  bar-and- 
beach  state  (Fig. 23c);  the  modal  state  appears  to  be  transitional  between 
the  two.  At  the  beginning  of  the  experiment  period,  the  beach  was  in  the 
longshore  bar-trough  state;  however  it  quickly  developed  into  the  rhythmic 
bar  and  beach  state  and  remained  in  that  state  for  the  duration  of  the 
period  even  though  the  trough  narrowed. 

Universally,  the  longshore  bar-trough  state  is  distinguished  by  its 
high  bar-trough  relief  combined  with  a  relatively  steep  beach  face.  Waves 
initially  break  over  the  bar  by  plunging  or  spilling  and  the  bar  region  is 
moderately  to  highly  dissipative  (i.e.  £  is  moderate  to  high).  In  contrast 
to  the  fully  dissipative  extreme,  however,  the  broken  waves  cease  their 
decay  after  passing  over  the  steep  inner  edge  of  the  bar  and  reform  within 
the  deep  trough.  The  much  steeper  beach  face  is  typically  reflective 
with  respect  to  the  partially  dissipated  and  subsequently  reformed  waves  and 
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c  is  commonly  in  the  vicinity  of  2.0.  Runup  is  consequently  relatively 
high.  In  terms  of  the  shore-normal  segregation  of  a  dissipative  bar  and 
reflective  beach  face,  the  rhythmic  bar  and  beach  state  is  similar  to 
the  bar  trough  state.  However,  rhythmic  longshore  undulations  of  the 
crescentic  bar  and  the  subaerial  beach  distinguish  the  rhythmic  bar 
and  beach  state  and  maintain  weak  to  moderate  rip  circulation,  the  rips 
are  persistently  located  in  embaymsits.  Spacings  between  shoreline 
protrusions,  crescentic  bar  horns  and  rips  are  typically  on  the  order  of 
200-300  metres  (Short,  1979  a  and  b,  Wright  et  al  1979  a  and  b; 

Wright  and  Short,  in  press)  . 

It  is  interesting  to  note  here  that,  despite  the  pronounced  temporal 
variability  of  deepwater  wave  conditions, a  characteristic  feature  of  the 
Bass  Strait  environment ,  temporal  variability  of  the  beach  state  is 
comparatively  small.  The  beach  and  surf  zone  appear  to  alternate  largely 
between  two  similar  states.  This  contrasts  markedly  with  the  tendency  for 
beaches  and  surf  zone  in  New  South  Wales  which  experience  dramatic  temporal 
variations  which  take  them  through  the  full  range  of  all  six  beach  states 
shown  in  Figure  23  (Short  1979  an  and  b;  Short  and  Wright  1981;  Wright 
et  al  1979  a  and  b;  Wright  and  Short,  in  press) .  The  comparative 
constancy  of  beach  state  (not  necessarily  of  shoreline  position)  which 
characterizes  the  Gippsland  beaches  is  probably  due  to  the  effects  of 
nearshore  frictional  dissipation  on  reducing  the  variability  of  breaker 
height  as  discussed  in  the  previous  section. 

Figure  24  shows  the  profile  configuration  of  the  beach  and  surf  zone 
of  Eastern  Beach  as  it  appeared  in  the  early  phases  of  the  experiment 
period.  Figure  25  a  and  b  show  the  plan  configurations  near  the  beginning 
and  end  of  the  experiment.  Aerial  oblique  and  ground  photographs  of  the 
study  site  are  shown  in  Figure  26  a  and  b.  Although  the  absolute  width  of 
the  trough  varied  considerably  over  the  experiment  period,  the  main  features 
as  shown  in  Figure  24  remained  consistent.  Waves  broke  over  the  shallow  bar 


by  plunging  or  spilling  and  breaker  heights  were  on  the  order  0.6  -1.2 


PRIMARY  BREAK 

PLUNGING  TO  SPILLING 
£,  =  30-100 


Figure  24.  Profile  and  major  dynamic  features  of  the  Eastern 
Beach  surf  zone.  Based  on  survey  of  5  May,  1981. 
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Figure  25.  Plan  configuration  of  the  bar,  trough  and  beach  at  the 
experiioent  site  and  locations  of  survey  lines 
(a)  5  May  1981  (b)  16  May  1981. 


metres.  The  inner  edge  of  the  bar  descended  abruptly  into  the  trough 

which  was  1.5  -  2.0  metres  deep.  The  beach  face  was  steep  and  maintained 

a  gradient  of  tan  P  =  0.10.  Over  the  step  at  the  base  of  the  beach 

there  was  a  violent  shore  break  (Fig. 27a)  followed  by  an  explosive 

surge  of  high  runup  on  the  beach  face  (Fig. 27b).  To  distinguish  the 

characteristics  of  the  outer  primary  break  and  the  shore  break  it  is 

necessary  to  define  two  scaling  parameters,  r  (eq.66) ,  an  outer  (or 

surf)  value  e  computed  from  the  primary  breaker  amplitude  and  the  local 
s 

qradient  of  the  bar  crest  and  a  beach  face  value,  r.,  ,  computed  on  the 

b 

basis  of  the  shore  break  amplitude  and  the  beach  gradient.  These 

parameters  had  values  of  f.  =  30-100  and  e,  =  1. 5-3.0  indicating,  dissi- 
*  s  b 

Dative  and  reflective  conditions  respectively  (Fin. 24). 

Consistent  with  observations  of  bar-trough  topoqraphies  elsewhere 
(e.g.  Moruya  Beach,  N.S.W.:  Wright  1982)  the  profile  between  the  upper 
limit  of  the  beach  face  and  the  inner  edge  of  the  bar  conformed  well  to 
an  exponential  model  of  the  form 

h  =  h  (1  -  e~ax)  0-7) 

00 

where  h  is  local  depth,  x  is  distance  seaward  and  h^  is  the  asymptotic 

depth  at  which  dh/dx  ~  0.  For  the  Eastern  Beach  case  shown  in  Figure  74 

h  =  1.7m  and  a  =  0.06.  The  distance  from  the  limit  of  runup  to  the 
00 

inner  edge  of  the  bar  (trough  width)  was  about  100  metres  during  the 
first  few  days  of  the  experiment  but  decreased  to  50  metres  as  the  bar 
migrated  toward  the  shore. 

In  plan  view,  the  beach  and  surf  zone  exhibited  long  crescentic 
bar  and  beach  rhythms  with  an  average  longshore  wave  lenqth  of  250  metres 
separating  shoreline  protrusions,  bar  crescents  and  rips  (Figs. 25  and 
26).  Rips  occupied  the  regions  of  shoreline  embayments  which  were  also 
the  regions  of  maximum  trough  width.  In  addition  to  the  large-scale 


longshore  rhythms,  the  steep  beach  face  exhibited  beach  cusps  spaced 
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at  intervals  of  about.  Ju  -  40  metres.  These  cusps  were  best  developed 
along  sections  fronted  by  the  most  pronounced  bar-trouqh  relief  and 
were  noorlv  developed  in  rip  bavs. 

''bsrrvcd  process  signature  of  the  surf  zone 

Table  6  summarizes  some  of  the  more  important  modes  of  fluid  motion 
which  contribute  to  the  shear  stresses  t  and  net  water  transport,  V  in 
surf  zones  in  qeneral.  These  modes  of  motion  include:  (1)  oscillatory 
flows  correspondinq  directly  to  the  incident  waves:  (2)  oscillatory  or 

quasi -oscillatory  flows  correspondinq  to  standing  waves  and  edge  waves 
at  frequencies  lower  than  incident  wave  frequencies,  particularly  those 
at  subharmonic  (period  twice  incident  wave  period)  and  infraqravity 
(period  >  30  secs)  frequencies:  and  (3)  net  circulations  Generated  by 
wave  enerqy  dissipation,  specifically  longshore  currents,  rips,  and  net 
return  flows.  Our  observations  of  different  surf  zones  show  that  the 
relative  velocity  magnitudes  of  these  different  modes  of  motion  depend  on 
morphodvnamic  state  and  also  vary  with  location  within  the  system 
(Wriatb  1982-  Wriqht  and  Short,  in  press). 

Tn  highly  dissipative  surf  zones,  shoreward  decay  of  incident 
waves  is  consistently  accompanied  by  shoreward  growth  of  infraoravi ty 
standinq  waves.  In  the  inner  surf  zone  of  such  beaches,  currents 
associated  with  infragravity  standing  waves  are  2  to  3  times  stroncier 
than  the  orbital  velocities  of  incident  waves  and  are  also  usually 
dominant  over  mean  currents  (Wright  et  al ,  1982).  On  reflective  beaches, incidc 
waves  and  subharmonic  edcre  waves  are  stronqly  dominant  and  mean  currents 
and  infraqravity  oscillations  are  very  weak  (Wriqht  et  al  1979a:  Wrioht 
1982:  Wriqht  and  Short,  in  press).  For  any  oiven  set  of  wave  conditions, 

the  strongest  rips  and  associated  feeder  currents  occur  in  association 
with  the  intermediate  transverse  bar  and  rip  state.  Tt  is  in  the  presence 
of  the  lonashore-bar- trough  and  rhvthmic  bar  and  beach  states  that  the 
most  complex  process  sionatu^es  prevail  (Wriqht  and  Short,  in  press: 


75 


Wright  1982).  Under  these  morphodynamic  conditions,  important  roles  are 
played  by  all  of  the  inodes  of  motion  listed  in  Table  6. 

In  a  succession  of  experimental  runs,  water  surface,  n  (pressure) 
and  horizontal  currents  u  and  v  time  series  were  recorded  at  different 
longshore  and  shore  normal  positions  within  the  surf  zone.  From 
statistical  and  spectral  analyses  of  these  time  series,  we  have  estimated 
the  velocity  amplitudes  corresponding  to  the  incident  waves  and  to  sub- 
harmonic  and  infragravity  standing  waves,  the  velocities  of  mean  currents, 
and  the  significant  heights  of  water  surface  oscillations  at  the  different 
frequencies.  The  results  of  these  analyses  from  several  of  the  more 
important  runs  are  listed  in  Table  7  together  with  the  positions  and 
depths  of  the  sensors.  From  an  amalgamation  of  the  data  in  Table  7  the 
process  signature  diagrams  shown  in  Figures  28  and  29  have  been  constructed. 

The  example  of  a  "process  signature  cross  section"  shown  in  Figure  28 

is  based  on  data  from  10  May  across  profile  6  after  the  bar  had  migrated 

shoreward.  The  upper  histograms  show  absolute  significant  heights  of 

oscillations  at  incident  wave  (Hg) ,  subharmonic  (Hsuij)  and  infragravity 

(H.  )  frequencies.  The  lower  histograms  show  the  corresponding  u  and  v 
in 

velocity  amplitudes.  It  is  apparent  from  Figure  28  that  motions  at  incident 
wave  frequency  are  dominant  at  all  locations.  However,  there  is  a  shoreward 
growth  of  subharmonic  and  infragravity  motions. 

The  composite  plan  view  process  signature  diagram  shown  in  Figure 
29  is  based  on  results  of  several  runs.  In  this  case,  velocity  amplitudes 
are  shown  on  dimensionless  scales  expressed  relative  to  the  magnitudes  of 
the  bottom  orbital  velocities  of  incident  waves.  The  general  surf  zone 
circulation  pattern  and  the  average  magnitudes  and  directions  of  time- 
averaged  currents  V  are  also  indicated.  Again,  the  incident  waves  are 
dominant  throughout  the  surf  zone  but  subharmonic  and  infragravity 
oscillations  become  pronounced  in  certain  regions.  Specifically  sub¬ 
harmonic  oscillations  increase  in  amplitude  with  proximity  to  the  beach 
within  the  trough  whereas  the  effects  of  the  infragrav  »  motions  are 
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Figure  28.  Process  signature  cross  section. 
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most  pronounced  in  the  rip  region.  We  will  discuss  each  of  the  processes 
separately . 

Incident  waves  and  breakers 

From  Table  7  two  important  points  concerning  the  behaviour  of  surf 

zone  oscillations  at  incident  wave  frequency  can  be  made.  The  first 

feature  deserving  mention  is  the  comparative  constancy  of  wave  heights 

from  the  inner  edge  of  the  bar,  across  the  trough  to  the  shore  break. 

(Fig.. 76b).  This  tendency  for  the  waves  to  maintain  their  height  after 

breaking  and  partially  dissipating  over  the  outer  portion  of  the  bar 

contrasts  markedly  with  patterns  observed  in  fully  dissipative  surf  zones 

(e.g.  Wright  et  al  1982).  In  low  gradient  dissipative  surf  zones,  surf 

height  decavs  progressively  shoreward  maintaining  a  constant  height-to- 

depth  (y  =  H,  /h)  ratio  (Wright  et  al  1962) .  However,  it  is  clear  from 
b 

Table  7  that  after  passing  over  the  bar,  the  surf  at  Eastern  Beach 
ceased  to  decay  and,  instead,  experienced  some  subsequent  shoaling  over 
the  inner  margins  of  the  trough  which  resulted  in  the  shore  break  being 
somewhat  higher  than  the  waves  over  the  inner  edge  of  the  bar  (Fig. 27a). 

The  second  point  is  that  the  height  of  the  shore  break  varied  by  only 
a  few  centimetres  throughout  the  study  period  despite  appreciable  temporal 
variations  in  deepwater  height,  H^.  Some  of  this  reduction  in  shore  break 
variation  was  undoubtedly  due  to  the  process  of  nearshore  frictional 
"filtering''  previously  described.  However,  much  of  the  filtering  was 
almost  certainly  due  to  the  process  of  partial  dissipation  over  the  bar. 

It  would  appear  that  this  pronounced  bar-trough  topography  is  effective 
in  limiting  the  heights  of  the  waves  which  actually  reach  the  beach  face. 

Standing  waves  and  edge  waves 

The  whole  or  partial  reflection  of  waves  from  the  shore  produces 
standing  waves  consisting  of  the  combination  of  incident  and  reflected 
waves.  Depending  in  part  on  the  degree  to  which  the  reflected  wave  is 
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refracted  as  it  propagates  into  deep  water,  standing  waves  may  be  either 
leaky  mode  standing  waves  (e.g.  Suhayda,  1974)  or  trapped  edge  waves 
(e.g.  Eckart,  1951;  Guta  and  Davis,  1974).  In  the  leaky  mode  case, 
reflected  energy  is  reradiated  back  to  sea  without  being  trapped  inshore 
and,  for  normally  incident  oscillations,  there  is  no  shore-parallel 
motion  or  amplitude  variation.  Alternatively,  reflected  energy  may  be 
trapped  inshore  by  refraction  to  produce  edge  waves  which  may  grow 
through  a  weak  resonant  instability  by  extracting  energy  from  incident 
waves  (Guza  and  Davis,  1974). 

The  theory  of  edge  waves  is  discussed  with  specific  reference  to 

beaches  by  Bowen  and  Inman  (1969);  Guza  and  Davis  (1974);  Ursell  (1951); 

and  Guza  and  Bowen  (1974) .  Edge  wave  radian  frequency,  w  (=2t/T 

e  e 

where  T  is  edge  wave  period)  and  alonashore  wave  number,  k  (=27r/L 
e  e  e 

where  L  is  alonashore  edge  wave  length)  on  a  beach  of  constant  slope  3 
e 

are  related  by  Eckart' s  (1951)  dispersion  relationship, 

a;  2  =  ok  (2n  +  1)  tan  6  (68) 

e  "  e 

where  n  is  edge-wave  mode  number  (n  =  0,1, 2, 3  ...),  designating  the  number 
of  nodal  positions  (zero  crossings)  in  the  x  direction  (perpendicular 
to  the  shore) .  For  the  case  of  an  exponential  inshore/beach  profile. 

Ball  (1967)  gives  the  edge  wave  dispersion  relationship  in  the  form, 

ui  2  =  4  qah  {  (2n  +  1)  (a2  +  4k  2)1*  -  a(2n2  +  2n  +1)}  (69) 

e  -  oo  e 

where  a  and  h  are  defined  by  equation  67. 

Normal  to  the  shore,  the  amplitudes  of  both  leaky  mode  standing 
waves  and  edge  waves  vary  from  a  maximum  at  the  shore  where  the  first 
antinode  occurs  to  zero  at  a  succession  of  nodes  which  are  separated 
by  antinodes  of  progressivelv  diminishing  amplitude.  The  variations  in 
amplitude  normal  to  the  shore  follow  the  offshore  dependence  functions 


which  are 
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(a)  Leaky  modes 

Ft  (x)  =  J  [2(  ■■<l,2-a  x/H 

L  o  [_  g  tan  S  J 


(Suhayda,  1974)  where  subscript  L  refers  to  leaky  mode  and  J  is  the 

o 

zero-order  Bessel  function 


(b)  Edge  waves 

F  (x)  =  exp  (-k  x)  Ln  (2k  x)  (71) 

e  e  e 

where  Ln  is  the  Laguerre  Polynomial  of  order  n  which  is  the  same  as  the 
edge  wave  mode  number. 

For  leaky  mode  standing  waves,  standing  edge  waves,  and  longshore- 
progressive  edge  waves  the  space  and  time  variations  of  surface  elevation 
n,  and  currents  u  and  v  as  expressed  by  the  single-valued  velocity 
potential  ♦(n--|||,u--^,v--|i) 

(a)  Leaky  mode 

a  =  aLg  Ft  (x)  sin  tot  (72) 

Li  -  Li  Li 

“l 

(b)  Edge  wave  standing  in  longshore  direction 

<t>  *  3eg  F  (x)  cos  k  y  cos  to  t  (73) 

e  -  e  e  e 

u 

e 

or 

(c)  Edge  wave  progressive  alongshore 

<P  =  aeg  F  (x)  sin  (k  y  -  u  t)  (74) 

e  -  e  e  e 

to 

e 

where  a^  and  a^  are  the  standing  wave  amplitudes  at  the  shoreline. 


I 
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The  simplest  edge  wave  model  is  that  of  the  synchronous  edge  wave  with 

frequency  equal  to  the  frequency  ok  of  the  forcing  incident  wind  waves 

or  swell  (Bowen  and  Inman,  1969;  Guza  and  Inman,  1975) .  The  synchronous 

edge  wave  is  also  the  highest  freauency  edge  wave  likely  to  occur  on  a 

beach.  Although  synchronous  edge  waves  have  been  produced  in  the  laboratory, 

they  have  not  been  shown  to  be  important  on  natural  beaches.  Guza  and 

Davis  (1974)  have  shown  that  it  is  more  common  for  edge  waves  with 

frequencies  lower  than  to  grow  by  a  weak  resonant  interaction  between 

the  incident  wave  and  two  edge  waves  of  frequencies  u;  and  w 

e(l)  e(2) 

travelling  in  opposite  directions  alongshore.  In  this  resonant  triad, 

the  conditions  m  +  <d  =  w,  and  L  ...  =  L  must  prevail. 

e(l)  e  (2)  l  e(l)  e(2) 

Guza  and  Davis  (1974)  and  Guza  and  Bowen  (1977)  found  that  the  most 

easily  excited  edge  wave  pair  is  that  of  subharmonic  edge  waves  of 

modes  n  =  0,0,  in  which  u  =  in  ...  =  w./2.  In  combination,  these 

e (1)  e (2)  l 

two  edge  waves  yield  a  standing  subharmonic  edge  wave  with  a  period  equal 
to  twice  the  incident  wave  period.  Our  observations  of  reflective 
beaches  in  Australia  indicate  that  subharmonic  edge  waves  are  normally 
well  developed  on  those  beaches  and  that  they  produce  cusps  spaced  at 
one-half  the  edge  wave  length  (Wright  et  al  1979;  Wright  1982;  Wright 
and  Short  in  press) . 

Subharmonic  oscillations  do  not  appear  on  highly  dissipative  beaches 
but  large-amplitude  infragravity  standing  waves  with  periods  of  100-200 
seconds  are  invariably  present  and  dominate  the  inner  surf  zones  of  those 
beaches  (Wright  et  al  1982;  Wright  1982) .  Recent  studies  by  Holman 
(1981)  and  Huntley  et  al  (1981)  suggest  that  these  infragravity  or 
"surf  beat"  oscillations  are  probably  edge  waves.  Observations  from 
intermediate  topographies  such  as  that  which  prevailed  at  Eastern  Beach, 
typically  reveal  lower  amplitude  and  higher  frequency  (period  =  30-60  secs) 
infragravity  standing  waves  as  well  as  subharmonic  oscillations  near  the 
beach  (Wright  et  al  1979;  Chappell  and  Wright  1979;  Wright  1982). 
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Power  spectra  and  cross  spectra  of  water  surface  elevation,  n, 
and  shore  normal  current,  u,  from  the  bar  trough,  and  step  regions  of 
the  surf-zone  are  shown  in  Figures  30-32.  Frequencies  which  are  standing 
as  interpreted  from  phase  angles  of  ir/2  between  n  and  u  are  indicated  by 
the  horizontal  bars.  The  spectra  show  dominance  of  energy  at  swell 
frequency  but  with  appreciable  standing  wave  energy  over  a  band  encom¬ 
passing  both  subharmonic  and  infragravity  frequencies.  From  Table  7  the 
subharmonic  oscillations  typically  had  heights  of  0.10  m  to  0.25  m  and 
associated  velocity  amplitudes  of  0.10  m  sec  ^  up  to  0.60  m  sec 
Infragravity  oscillations  had  heights  of  0.15  m  to  0.40  m  and  velocity 
amplitudes  (u  and  v)  on  the  order  of  0.20  m  sec  1  to  0.35  m  sec  \ 
However,  as  would  be  expected  for  standing  waves  and  for  edge  waves 
standing  alongshore,  and  hence  exhibiting  alternate  nodes  and  antinodes, 
the  heights  and  velocity  amplitudes  of  the  infragravity  and  subharmonic 
oscillations  varied  considerably  with  location  in  the  surf  zone  as  shown 
in  Figures  28  and  29. 

The  subharmonic  oscillations  were  identified  as  such  on  the  basis 
of  their  having  significant  standing  wave  energy  at  frequencies  close 
to  one  half  the  frequency  of  the  incident  swell.  These  frequencies  were 
typically  in  the  range  of  0.038  -  0.056  hz  (periods  of  18  to  26  seconds) . 
It  must  be  noted,  however,  that  there  was  not  always  a  clear  separation 
between  what  we  refer  to  as  subharmonic  and  the  longer  period  swell  or 
between  "subharmonic"  frequencies  and  the  higher  frequency  infragravity 
motions.  There  was  in  fact,  a  broad  band  of  standing  wave  energy  from 
0.056  hz,  through  to  about  0.015  hz.  Having  expressed  this  caveat,  we 
will  proceed  on  the  assumption  that  the  subharmonic  oscillations  are 
edqe  waves  since  observations  from  other  reflective  and  intermediate 
beaches  suggest  that  they  are  (Wright  et  al  1979?  Wright  1982) . 

From  Figure  29  the  velocity  amplitudes  associated  with  subharmonic 
frequencies  are  seen  to  decay  seaward  from  a  maximum  above  the  step 
indicating  that  these  oscillations  are  largely  confined  to  the  vicinity 
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of  the  beach  face.  This  is  roughly  consistent  with  expectations  for  r<_  le 
zero  edge  wave  which  is  the  most  probable  for  a  reflective  beach  face 
(Guzaand  Davis,  1974).  For  a  24  second  (0.042  hz)  subharmonic  n  =  0  edge 
wave  on  the  observed  beach  face  slope  of  0.10,  equation  68  gives  an  edge 
wave  length  of  90  metres.  This  is  close  to,  but  exceeds,  twice  the 
observed  cusp  spacing  of  35  to  40  metres.  An  edge  wave  with  a  period  of 
22  seconds  would  have  an  alongshore  length  of  76  metres  which  is  very 
close  to  twice  the  cusp  spacing.  However,  Wright  (1982)  has  found  that  on 
reflective  and  bar-trough  beaches  where  the  beach  face  and  inner  surf 
zone  profile  fits  the  exponential  model  expressed  in  equation  67  the 
cusps  spacing  fits  most  closely  to  the  edge  wave  lengths  predicted  on  the 
basis  of  Ball's  (1967)  equation  (eq.69) .  On  this  basis  using  the  observed 
a  ~  0.06  an'  h  =  1.7  m  a  24  second,  n  =  0  edge  wave  would  have  a  length  of 

CO 

73  metres.  Regardless  of  whether  we  use  equation  68  or  69  it  appears  very 
likely  that  the  cusps  which  occurred  on  the  reflective  beach  face  in  the 
lee  of  the  offshore  bar  were  formed  by  subharmonic  edge  waves  standing 
alongshore.  This  is  consistent  with  numerous  observations  on  natural 
reflective  and  intermediate  beaches  (Huntley  and  Bowen,  1979;  Wright 
et  al  1979;  Wright  1982)  and  in  the  laboratory  (Guza  and  Inman  1975)  that 
standing  subharmonic  edge  waves  form  cusps  spaced  at  one  half  the  edge 
wave  length  when  conditions  are  sufficiently  reflective. 

The  infragravity  standing  waves  in  the  Eastern  Beach  surf  zone 
occurred  at  much  higher  frequencies  and  were  of  considerably  lower 
amplitude  than  those  commonly  observed  in  highly  dissipative  surf  zones 
(e.g.  Wright  et  al  1982).  Typical  periods  of  the  infragravity  oscillations 
at  Eastern  Beach  were  on  the  order  of  35-60  seconds  (0.028  -  0.017  hz) . 

These  frequencies  are  very  similar  to  those  observed  in  the  presence  of 
bar-trough  topographies  elsewhere  (e.g.  Moruya  Beach  N.S.W.:  Wright  et  al 
1979;  Chappell  and  Wright  1979;  Wright  1982) .  These  intermediate 
frequency  infragravity  standing  waves  have  been  discussed  by  Chappell  and 
Wright  (1979)  and  Wright  (1982)  who  suggest  that  the  frequencies  may  be 
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constrained  or  selected  by  the  dimensions  of  the  trough.  This  suggestion 
has  some  credibility  in  the  Eastern  Beach  case  since  the  natural  frequency 
of  the  trough  is  estimated  to  be  on  the  same  order  as  the  observed  infra¬ 
gravity  frequencies. 

It  is  tempting  to  suggest  that  the  rhythmic  spacings  of  the  larger 
scale  surf  zone  features  such  as  rips,  crescentic  bar  "horns"  and  shoreline 
protrusions  are  related  to  long  infragravity  edge  waves.  We  presented 
arguments  for  this  in  earlier  papers  (Wright  et  al  1979;  Wright  1982) . 
Recent  work  by  Holman  and  Bowen  (1982)  has  shown  how  rhythmic  crescentic 
bar  topography  such  as  that  which  prevails  at  Eastern  Beach  as  well  as  more 
complex  quasi-rhythmic  surf  zone  topographies  can  be  explained  by  drift 
patterns  related  to  the  interactions  between  two  edge  wave  modes .  The 
simple  crescentic  bar  and  rhythmic  beach  pattern  like  that  at  Eastern 
Beach  is  considered  to  represent  an  equilibrium  response  to  a  standing 
edge  wave  composed  of  two  edge  wave  modes  of  the  same  frequency  travelling 
in  opposite  directions  (Holman  and  Bowen,  1982) .  This  will  produce 
longshore  topographic  periodicities  with  wave  lengths  of  one  half  the 
length  of  the  individual  component  edge  waves.  For  the  observed  averaged 
surf  zone  beach  gradient  of  about  0.03,  a  60  second,  mode  1  edge  wave 
in  the  Eastern  Beach  surf  zone  should,  on  the  basis  of  equation  68,  have 
had  an  alongshore  wave  length  of  about  500  metres  which  is  close  to  twice 
the  alongshore  wave  length  of  the  larger  scale  surf  zone  features.  Such 
edge  waves  must  extend  well  seaward  of  the  bar  crest. 

A  comparison  of  pressure  (or  n)  spectra  from  the  step,  trough  and 
bar  regions  of  the  surf  zone  is  shown  in  Figure  33.  The  important  features 
are  the  sharp  infragravity  peaks  over  the  step  and  bar  accompanied  by  an 
energy  trough  at  the  same  frequency  in  the  centre  of  the  topographic  trough. 
This  suggests  the  presence  of  shore  normal  antinodes  at  the  beach  and  over 
the  bar  with  the  intervening  node  located  in  the  trough.  This  is  consistent 
with  the  logical  and  widely  held  notion  that  troughs  should  be  best  formed 
in  the  presence  of  the  stronger  horizontal  currents  of  the  nodal  region 


SPECTRAL  DENSITY 
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Figure  33.  Power  spectra  of  water  surface  elevation,  n,  time  series 
recorded  simultaneously  over  the  bar,  in  the  trough  and 
over  the  step,  At  =*  1.5  sec,  df  -  12.  Note  that  infra¬ 
gravity  peaks  over  the  bar  and  over  the  step  are 
acconpanied  by  an  energy  "trough"  at  the  same  frequency 
within  the  topographic  trough. 
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whereas  bars  should  accumulate  beneath  antinodes  where  drifts  are  weaker. 
In  the  shore  parallel  dimension,  both  u  and  v  components  of  the  infra¬ 
gravity  motions  intensify  in  the  embayment  region  where  the  rip  is  located 
(Fig. 29;  Table  7) ,  suggesting  that  this  is  probably  the  location  of  the 
shore  parallel  node.  A  comparison  of  the  possible  instantaneous  water 
surface  topographies  associated  with  the  co-existing  subharmonic  and 
infragravity  edge  waves  is  shown  in  Figure  34. 

Surf-zone  circulation  and  net  currents 

Hie  net  (time  averaged)  circulation  features  observed  in  the  Eastern 
Beach  surfzone  were  of  three  general  types  (1)  vertically  segregated  shore 
normal  drifts;  (2)  horizontally  segregated  rip  circulation;  and 
(3)  longshore  currents.  The  vertically  segregated  drifts  were  prominent 
over  the  bar  and  within  the  trough  in  regions  away  from  rips.  This  form 
of  circulation  was  similar  to  that  described  for  the  dissipative  surf  zone 
at  Goolwa,  South  Australia  by  Wright  et  al  (1982)  and  was  characterized  by 
net  onshore  transport  in  the  surface  layer  and  net  offshore  transport 
near  the  bottom.  Onshore  wind  stress  probably  played  in  a  role  in  main¬ 
taining  this  form  of  circulation.  '.»e  nature  of  this  circulation  is 
illustrated  by  the  results  of  a  vertical  profile  of  shore-normal  (u) 
velocity  measured  in  the  trough  on  profile  6  at  a  depth  of  2.0  metres. 
Figure  35  shows  the  time  series  of  u,  filtered  with  a  100  second  filter 
to  remove  wave  effects,  from  0.04  m  0.22  m  0.45  m  and  0.91  m  above  the 
bed.  Closest  to  the  bed,  flow  is  seen  to  be  exclusively  offshore  whereas 
it  is  predominantly  onshore  on  0.91  m  above  the  bed.  Onshore  transport 
prevailed  throughout  the  upper  metre  of  the  water  column.  As  indicated 
by  the  values  of  time  averaged  near  bottom  velocities  (<u>)  in  Table  7, 
flows  below  0.30  m  above  the  bed  were  seaward  throughout  the  study  period. 

Since  onshore  transports  in  the  upper  half  (or  more)  of  the  water 
column  were  substantially  stronger  than  seaward  transports  near  the  bed, 
it  may  be  concluded  that  the  depth  averaged  transport  was  shoreward  in 
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Figure  35.  Shore-normal  current  (u)  time  series  from  4  elevations 
showing  vertical  segregation  of  onshore  (upper  level) 
and  offshore  (lower  levels)  flows,  an  electronic  filte. 
with  a  100  second  time  constant  has  been  applied  to 


remove  incident  waves. 
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the  vicinity  of  profile  6  which  was  situated  near  a  shoreline  protrusion. 
This  excess  onshore  transport  was  accommodated  by  the  net  horizontal 
circulation  pattern  shown  in  Figure  29.  Water  flowed  alongshore  within 
the  trough  and  turned  offshore  as  a  wide  rip  in  the  region  of  profile  4. 
Within  the  rip,  flow  was  seaward  at  all  depths.  Unfortunately,  we  were 
unable  to  obtain  a  vertical  velocity  profile  owing  to  the  strong  current 
drag  on  the  vertical  mast.  However,  near  bottom  (0.20m  above  the  bed) 
velocities  of  up  to  0.30m  sec  ^  were  recorded  and  dye  tracking  of  surface 
flows  indicated  stronger  surface  flows.  However,  the  observed  rip 
velocities  were  substantially  weaker  than  those  which  normally  occur  in 
association  with  the  transverse  bar  and  rip  type  of  topography  (Wright 
and  Short  in  press) .  Despite  the  fact  that  the  rip  was  not  excessively 
strong,  it  was  nevertheless  effective  in  causing  seaward  bedload  transport 
of  sand  as  evidenced  by  the  presence  of  seaward-migrating  dune  bedforms. 

For  a  few  days  during  the  experiment  period,  strong  longshore  currents 
prevailed  within  the  trough  related  in  part  to  oblique  wave  incidence  and 
in  part  to  the  "feeding"  of  the  downdrift  rip.  We  succeeded  in  obtaining 
vertical  velocity  profiles  of  these  currents,  some  examples  of  which  are 
shown  in  Figure  36.  Although  longshore  currents  over  simple  topography 
have  been  the  subject  of  considerable  attention  (e.g.  Komar  1975;  1976; 

Longuet-Higgins,  1970a  and  b)  there  is  very  little  field  data  on  the  nature 
of  longshore  currents  within  pronounced  troughs  or  on  the  vertical  current 
structure.  Some  more  detailed  analysis  and  discussion  of  our  data  is 
therefore  warranted. 

Time  averaged  longshore  current  velocities  (<v>)  near  the  surface 
exceeded  0.5m  sec  1  in  some  runs  while  as  close  as  0.075m  above  the  bed 
(above  the  level  of  ripple  troughs)  velocities  were  still  as  high  as 
0.25m  sec  ^  (Fig. 36).  Although  these  longshore  velocities  were  on  the 
same  order  of  magnitude  as  the  wave  orbital  velocities,  the  bottom  boundary 
layer  was  completely  dominated  by  the  waves  which  acted  in  the  shore 
normal  direction  orthogonal  to  the  longshore  currents.  This  is  indicated 
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LONGSHORE  VELOCITY  PROFILE  .  MAY.  1982 
EASTERN  BEACH.LAKES  ENTRANCE. G/PPSLAND. 

A  10  SEC  FILTER, SAMPLED  EVERY  6  SECS. 

«  100  SEC  FILTER, SAMPLE  INTERVAL  1  MIN. 

□  100  SEC  FILTER, SAMPLE  INTERVAL  1  MIN. 
x  FAST, SAMPLE  INTERVAL  6  SECS. 


Figure  36.  Vertical  profile  of  time-averaged  longshore  current 
velocity,  <v>,  as  recorded  within  the  trough  of 
Eastern  Beach  in  a  succession  of  runs.  The  position  of 
the  lowest  flow  meter  is  shown. 
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by  two  lines  of  evidence:  (1)  the  bed  forms  were  fully  developed, 
sharp-crested  wave  induced  ripples  (heights,  nr  =  0.08m)  with  long 
crests  aligned  shore  parallel;  and  (2)  vertical  profiles  of  suspended 
sediment  concentration  (discussed  further  in  the  next  subsection)  were 
perfectly  exponential  similar  to  those  observed  under  pure  wave  motion 
and  concentration  magnitudes  were  not  significantly  larger  than  would 
be  expected  from  waves  alone . 

On  the  other  hand,  it  appears  that  the  longshore  current  profiles 
were  strongly  influenced  by  the  vortices  generated  by  the  wave  oscillations 
over  the  ripples  since  they  were  not  at  all  logarithmic  (Fig. 36).  Between 
the  elevations  of  z  =  0.07m  and  0.30m  above  the  bed,  which  is  the  domain  of 
the  vortices  shed  by  the  ripples,  the  longshore  current  velocity  was 
nearly  constant.  This  is  understandable  since  the  time  averaged  longshore 
shear  stresses,  t  must  be  staticly  determined,  as  in  steady  flows,  the 
velocity  gradients,  d  <v>/dz,  must  adjust  to  the  total  eddy  viscosity, 
vt<  that  is 


=  T/pv 


and  hence  be  small  where  is  large.  The  wave  induced  vortices  cause  a 

very  large  exchange  of  momentum  for  a  given  shear  rate,  and  thus  produce 

a  large  eddy  viscosity.  This  strong  effect  of  the  vortices  must  be 

related  to  the  fact  that  the  axes  of  the  vortices  are  oriented  parallel 

to  the  direction  of  the  longshore  current.  The  boundary  layer  structure 

is  apparently  entirely  different  in  rip  currents  where  wave  oscillations 

parallel  the  mean  flow  producing  a  "roller  effect".  In  order  to  model  the 

current  flow  near  the  bed,  we  will  consider  the  mean  longshore  flow  as  a 

steady  current  which  only  interacts  with  the  waves  by  way  of  the  eddy 

viscosity,  v ,  which  consists  of  two  components  (vt  =  vw  +  vc^ :  a  component, 

v  attributable  to  the  waves  and  a  component  v  attributable  to  the  current, 
w  c 


L 


Equation  75  relating  the  time  average  longshore  velocities  to  shear 


stress  then  becomes 


P 


(v 


w 


+  V  ) 

c 


d<v> 

dz 


=  T  (z) 


(76) 


The  eddy  viscosities,  v  and  v  are  related  to  vertical  velocity  fluctuations 

w  c 

only  and  are  thus  not  vectors  in  the  x-y  plane. 

We  assume  the  distribution  of  shear  stress  to  be  linear 


T(z)  =  t  (1-z/h) 
c 


(77) 


where  is  the  longshore  shear  stress  at  the  bed.  The  contribution  of 
the  longshore  current  to  eddy  viscosity  is 


v  =  C  v  z ( 1-z/h) 
c  c  *c 


where  v*c  is  the  longshore  friction  velocity 


(78) 


(79) 


For  "pure"  steady  flow,  the  constant  is  equal  to  von  Kantian's 
constant  (ic  =  0.4) . 

The  wave-induced  eddy  viscosity  can  be  evaluated  from  the 
empirical  formula 

_ CW  U*w  2 _  (80) 

w  *  l+1.34(*5f  )•  £exp(z/6) 
w  <5 

where  f^  is  the  friction  factor,  6  is  boundary  layer  thickness  and  u^ 
is  wave  friction  velocity.  Equation  80  was  obtained  by  Lundgren  (1972) 
from  analyses  of  boundary  layer  measurements  subsequently  published  by 
Jonsson  and  Carlsen  (1976). 
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We  note  that  is  initially  proportional  to  z  but  then  decays  exponentially 

with  z  farther  frcm  the  boundary  with  a  length  scale  of  decay  equivalent 

to  the  boundary  layer  thickness  6 . 

Exponential  decay  was  observed  by  Nakato  et  al  (1977)  and  MacDonald 

(1977)  in  measurements  of  vertical  velocity  fluctuations,  w  ,  in 

rms 

oscillatory  boundary  layers  (Fig. 37) .  Their  measurements  also  show 

that  the  vertical  length  scale,  l,  of  suspended  sediment  concentration 

“Z/Jl 

profiles  (c(z)  =  C  e  ;  where  C  is  concentration  at  the  bed)  are 
o  o 

identical  to  those  of  the  velocity  distribution.  The  correspondence 

between  w  and  c  is  shown  in  Figure  37.  Although  we  have  no  measurements 
rms 

of  w  ,  we  do  have  direct  field  measurements  of  c(z)  (discussed  further  in 
rms 

the  next  subsection)  obtained  simultaneously  with  the  longshore  current 
profiles.  We  can  use  these  data  to  obtain  an  estimate  of  the  length  scale, 
l,  (since  c(z)  showed  consistent  exponential  decay  with  z)  from 


l 


d  In  c  1 
dz 


0.093m 


(81) 


The  length  scale,  l,  can  then  be  used  to  replace  5  in  equation  80. 

By  replacing  6  with  Z,  and  inserting  and  as  obtained  from 
(78)  and  (80)  into  equation  76  we  get 


d<v> 

dz 


v?  (1-z/h) 
wc 


V  +  V 
w  c 


or 


<v2>  -  <v}> 


v^f  (1-z/h) 
"c 


V  +  V 

w  c 


dz 


(82) 


(83) 


where  <Vj>  and  <v2>  are  the  longshore  velocities  at  two  elevations,  Zt 

and  zj  respectively.  If  and  Cc  are  known,  equation  82  will  determine 
v*c  from  any  pair  of  measured  velocities.  However,  we  cannot  consider 
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Vertical  distribution  of  vertical  velocity,  H^,  and 
suspended  sediment  concentration  C  after  Nakato  et  al 
(1977)  aid  MacDonald  (1977) .  This  graph  indicates  that 
c,  when  known,  provides  an  acceptable  surrogate  for  w 


Fioure  37. 
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C  =c  =  k  =  0.4  since  this  will  greatly  underestimate  the  relative 
w  c 

importance  of  v  and  lead  to  a  poor  prediction  of  the  velocity  profile 
w 

as  illustrated  in  Figure  38.  However,  if  we  set  C  =  0.4  and  adjust 

w 

C  and  v  a  good  fit  is  obtained  for  C  =0.1  and  v.  =  0.24m  sec  \ 
c  *c  c  *c 

Alternatively,  as  shown  in  Figure  38,  an  equally  good  fit  is  found  for 

the  combination  C  =  0.4,  C  =  7.4  and  v.  =  0.10m  sec  \  The  resulting 
c  w  *c 

profiles  of  v  ,  v  ,  and  v  are  shown  in  Figure  39. 
wet 

Our  available  empirical  basis  is  insufficient  to  permit  reliable 

prediction  of  longshore  current  profiles.  However,  the  qualitative 

inplications  are  that  wave  generated  vortices  shed  from  ripples  have 

a  paramount  influence  on  longshore  current  profiles.  The  longshore 

bed  shear  stresses  corresponding  to  the  two  alternative  parameter  sets 

(Cw,  C^,  and  v^)  described  above  differ  by  a  factor  of  -17.  Independent 

measurements  of  t  ,  for  example  by  direct  measurements  of  mean  water  surface 
c 

slope  (dn/dy)  could  readily  indicate  which  set  of  parameters  is  most 
appropriate. 

Sediment  suspension  in  the  surf  zone 

The  observed  presence  of  wave-formed  ripples  within  the  trough 

together  with  the  observed  tendency  for  wave  oscillations  to  dominate 

the  boundary  layer  suggests  that  oscillatory  flows  at  incident  wave 

frequency  should  dominate  the  entrainment  and  suspension  of  sediment. 

As  discussed  previously  in  relation  to  sediment  transport  seaward  of  the 

surf  zone,  recent  laboratory  and  field  studies  indicate  that  oscillatory 

flows  over  rippled  beds  cause  predominantly  suspended  load  transport 

(albeit  close  to  the  bed)  and  that  suspended  sediment  concentrations  can 

be  modelled  as  functions  of  bottom  orbital  velocities  (a  cu) ,  sediment 

s 

fall  velocity,  w  ,  and  ripple  height  n  (Nielsen  et  al  1978;  Nielsen 
s  r 

1979;  Neilsen  and  Green  in  press  ,  Nielsen  in  preparation). 
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Figure  38.  Predicted  and  observed  longshore  velocity  profiles 


Figure  39.  Vertical  distribution  of  total  Cvfc) ,  wave  induced  (v^) 
and  current  induced  (v  )  eddy  viscosity. 


The  conventional  method  of  expressing  suspended  sediment 


concentration  profiles  is  in  terms  of  the  one-dimensional  diffusion 
equation 
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dc 

dz 
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w  c 
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0  or 


d  In  c 
dz 


(84) 


where  E  is  diffusivity,  c  is  suspended  sediment  concentration  and  w^ 
is  sediment  fall  velocity  (e.g.  Nielsen  1979;  Nielsen  and  Green  in 
press) .  However,  recent  work  by  Nielsen  (in  preparation)  suggests  that 
a  more  meaningful  approach  is  to  consider  c  in  terms  of  the  vertical 
length  scale,  SL,  in  the  form 


C(z)  =  C  exp  ( -z/i ) 
o 


(85) 


where  C  is  the  concentration  at  the  ripple  crest  level  (z  =  0)  and 
o 

l  is  as  defined  by  equation  81  that  is 


f  d  In  cl 

■  L""^j 


and  is  related  to  the  diffusivity  by 


Hz) 


E(z) 

w 

s 


(87) 


Our  change  in  formulation,  that  is  from  E  to  £,  is  inspired  by 
observations  of  concentration  profiles  for  different  sand  sizes  in 
the  same  flow.  A  set  of  such  observations  is  shown  in  Figure  40.  We 
see  that  at  least  the  upper  parts  of  the  profiles  tend  to  be  parallel 
although  the  settling  velocities  vary  by  a  factor  of  7.  If  the  entrain¬ 
ment  process  was  diffusive  and  the  diffusion  equation  was  strictly  valid, 
then  the  profile  slopes  should  be  inversely  proportional  to  the  settling 
velocity.  This  is  clearly  not  the  case.  Instead,  there  is  a  strong 
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Figure  40.  Vertical  suspended  sediment  concentration  profiles  for 
different  grain  sizes  (from  direct  field  measurements) 
Note  that  despite  appreciable  variation  in  settling 
velocity,  w  ,  the  upper  parts  of  the  profiles  are 
parallel  suggesting  that  the  suspension  of  the  various 
sizes  is  related  to  a  single  length  scale,  l. 
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indication  that  the  flow,  in  fact,  defines  a  single  length  scale,  l, 
appropriate  to  all  grain  sizes.  As  pointed  out  previously  in  connection 
with  Figure  37,  this  sane  length  scale  also  applies  to  the  decay  of 
vertical  velocity  fluctuations  with  elevation. 

Nielsen  (1979)  and  Nielsen  and  Green  (in  press)  have  shown  that 
Cq  is  a  function  of  the  skin  friction  Shield's  parameter 


O'" 


t  /(p  -  p)  go 
s  s 


h_ 

(p 


pf  (a  (ii)2 

w  s 

s  "  p)  gD 


s 


(88) 


where,  as  previously  f^  is  the  friction  factor,  a^  is  orbital  semi- 

excursion,  w  is  wave  radian  frequency  p  and  p  are  sediment  and  water 

s 

densities  and  D  is  grain  diameter.  The  length  scale,  l,  relative  to 
s 

ripple  hei^it,  n  ,  is  a  function  of  (a^u/w^)  where  w^  is  sediment  fall 
velocity  (Nielsen  in  preparation). 

Using  the  suction  device  described  in  the  methods  section,  we 
obtained  12  concentration  profiles  from  the  Eastern  Beach  surf  zone. 
Examples  of  concentration  profiles  from  separate  runs  are  shown  in 
Figure  41.  The  results  of  these  measurements  are  strictly  comparable  to 
field  results  from  other  beaches  as  well  as  laboratory  results  (e.g. 
Nielsen  and  Green  in  press)  and  indicate  that  the  concentration  profiles 
can  be  reasonably  modelled  on  the  basis  of  suspension  by  waves  alone. 
Estimates  of  Cq  based  on  the  Eastern  Beach  data  together  with  other  sets 
of  field  data  indicate  that  Cq  is  well  predicted  from 

C  =  .00018  O'3  (89) 

o 

A  comparison  of  the  predictive  curve  of  equation  89  with  field  data  is 
shown  in  Figure  42.  Data  points  from  the  Eastern  Beach  surf  zone  are 
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indicated  and  are  seen  to  conform  well  to  the  general  model. 
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Figure  41.  Suspended  sediment  concentration  profiles  from  Eastern 
Beach. 
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Figure  42.  C  as  a  function  of  0' . 
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The  length  scale,  l,  relative  to  ripple  height,  n  #  gives  the 
best  fit  to  the  relationship 


=  1.43  -  1.25  exp  ^-.0011 

from  which  we  can  predict  l  if  we  know  or  can  estimate  ripple  height, 

n  .  Nielsen  (1981)  showed  that  n  is  largely  the  function  of  the 
r  r 

mobility  number. 
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Available  data  suggests 


a 
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_ 0.23 _ 

1+0.022  exp  (4716) 


(92) 


Figure  43  shows  the  relationship  between  values  of  l  as  predicted  from 
equations  90  and  92  and  values  observed  in  the  laboratory  and  in  the 
Eastern  Beach  surf  zone.  Although  there  is  appreciable  scatter,  the  fit 
is  reasonably  good.  We  conclude  that  despite  the  complexity  of  processes 
in  the  presence  of  the  longshore  bar-trough  and  rhythmic  bar-and-beach 
morphodynamic  states,  suspended  sediment  concentrations  can  probably  be 
modelled  in  terms  of  entrainment  by  incident  waves  alone.  More  data  are 
still  needed  to  substantiate  this. 

Unfortunately,  we  still  lack  confidence  in  our  ability  to  predict, 
quantitatively,  the  rate  of  net  transfer  of  sediment  from  one  region  of 
the  surf  zone  to  another.  To  a  degree,  at  least,  this  obviously  involves 
superimposition  of  the  net  flow  field  on  the  suspended  sediment  concentration 
field.  This  is  reasonable  to  assume  for  the  alongshore  transport  where 
net  flow  acts  perpendicular  to  wave  oscillations.  As  was  the  case  over 
the  nearshore  profile,  longshore  advection  of  suspended  load  should  be 
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However,  as  explained  in  the  discussion  of  nearshore  processes,  we 
cannot  extend  this  argument  to  the  shore-normal  case  where  wave 
oscillations  and  net  flows  act  parallel  to  each  other.  As  we  pointed 
out  earlier,  phase  dependency  between  velocity  maxima  and  concentration 
maxima  can,  in  fact,  result  in  sediment  moving  opposite  to  the  direction 
of  the  net  flow.  This  point  is  particularly  well  borne  out  by  the 
Eastern  Beach  case.  Although  net  near-bottom  flow  of  water  was  exclusively 
seaward  over  the  study  period  (e.g.  Fig. 35;  Table  7)  successive  surveys 
(discussed  in  the  following  section)  indicated  a  significant  net  onshore 
transport  of  sediment  over  the  same  period. 
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MORPHODYNAMIC  VARIABILITY 
OF  THE  BEACH  AND  SURF  ZONE 


Beach  and  surfzone  morphologies,  and  hence,  coupled  sets  of 
morphodynamic  processes,  represent  responses  to  prevailing  and 
antecedent  breaker  conditions  in  relation  to  the  sediment  composing 
the  morphology.  Variations  in  either  breaker  conditions  or  sediment 
size  can  be  expected  to  produce  variations  in  morphodynamic  state 
and  in  profile  configuration.  Discussion  of  the  relationships 
between  morphodynamic  state  and  wave  conditions  for  specific 
beaches  have  been  presented  by  Short  (1979  a  &  b;  1981)  ,  Wright  and 
Short  (in  press)  and  Wright  et  al  (1979) .  Goldsmith  et  al  (1982) 
have  considered  the  relationships  between  changing  bar-trough  and 
crescentic  bar  topographies  and  changing  wave  conditions.  We  are 
concerned  here  with  both  spatial  variability  along  the  coast  and 
temporal  variability  at  the  Eastern  Beach  site. 

Spatial  Variability 

Spatial  variability  in  morphodynamic  state  can  be  expected 
to  result  from  systematic  alongshore  variations  in  breaker  height 
(for  example  as  a  result  of  variations  in  nearshore  wave  modification 
or  in  exposure)  or  to  longshore  variations  in  sediment  size,  or  both. 
Spatial  variability  along  the  Gippsland  coast  was  assessed  from 
aerial  photos,  three  sets  of  aerial  reconnaissances, and  a  systematic 
set  of  surveys  and  sediment  sampling  at  35  sites  along  the  entire 
length  of  coast.  The  survey  and  sediment  sampling  programme  was 


conducted  in  December,  1981  over  a  short  period  during  which  deepwater 


wave  conditions  remained  constant,  eliminating  possible  temporal 
variations.  The  surveys  extended  from  the  foredune  to  the  outer 
bar.  The  locations  of  the  transects  are  shown  in  Figure  44 
and  the  resulting  profiles  are  shown  in  Figure  45.  Beach  states 
were  recorded  at  each  site. 

The  most  pronounced  spatial  variability  occurs  in  the  extreme 
northeastern  sector  between  Cape  Howe  and  Cape  Everard  (Figure  44) 
owing  to  the  combined  effects  of  variations  in  coastal  orientation 
and  sheltering  by  headlands.  Along  that  portion  of  the  coast, 
longshore  bar- trough  states  prevail  in  the  most  exposed  localities 
but  in  sheltered  areas  conditions  become  more  reflective.  However, 
from  Cape  Everard  southwest  to  Seaspray  (Profile  30)  beach  state 
remains  consistently  in  the  longshore-bar-trough  and  rhythmic-bar- 
and-beach  state  indicating  that  conditions  observed  at  the  Eastern 
Beach  site  are  typical  of  most  of  the  coast.  Between  Seaspray  and 
Reeves  Beach,  beach  state  begins  to  shift  toward  the  transverse-bar 
and  rip  state  as  energy  decreases.  South  of  Shallow  Inlet,  the 
sheltered  environment  favours  the  prevalence  of  reflective  states. 

Despite  the  comparative  constancy  alongshore  of  beach  state, 
particularly  between  Cape  Everard  and  Seaspray,  there  are  significant 
variations  in  the  profiles  as  shown  in  Figure  45  .  Although  the 
general  features  of  the  profile,  as  discussed  in  the  previous  section 
remain  consistent,  trough  width  and  the  slope  of  the  ir.’er  beach 
face  vary  considerably,  apparently  in  response  to  variations  in 
local  sediment  size.  In  general  the  beach  face  is  composed  of 
medium  sand  CtO.39  mm),  the  coarsest  material  (tl.  5  mm)  comprises 
the  step  at  the  base  of  the  beach  face ,  and  the  bar  is  composed  of 
finer  (,'o.32  mm)  sand.  However,  significant  coarsening  of  material 
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Figure  44.  Locations  of  beach  profiles  and  sediment  samples,  Deceirfcer,  1981. 
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Figure  45.  Spatial  variability  of  beach  profiles,  December,  1981 
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Figure  46.  Variation  of  beach  face  slope,  tan  8,  with  grain  size. 
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in  certain  areas,  notably  between  Stockyand  East  and  Seaspray  (Figure  44) 
causes  a  narrower  trough  and  a  steeper  subaerial  beach.  The  relation¬ 
ship  between  beach  face  slope  and  grain  size  is  particularly  obvious, 
consistent  with  the  observations  of  numerous  others.  This  relation¬ 
ship  is  shown  in  Figure  46. 

Temporal  Variability 

The  temporal  variability  of  a  beach  expresses  the  temporal 
variability  of  breaker  conditions.  On  the  New  South  Wales  coast, 
where  the  highly  variable  deepwater  wave  climate  experiences  only 
minor  nearshore  attenuation, breaker  conditions  are  highly  variable 
and,  accordingly,  exposed  beaches  experience  changes  which  can  take 
them  through  the  full  range  of  beach  states  (Figure  23  ;  Short 

1979;  1981;  Wright  et  al  1979  a  &  b) .  Although  the  deepwater  wave 
climate  of  Bass  Strait  is  at  least  as  variable  as  that  affecting 
the  New  South  Wales  coast,  we  have  shown  earlier  in  this  report  that 
nearshore  dissipation  of  waves  approaching  the  Gippsland  coast 
substantially  reduces  the  temporal  range  of  breaker  conditions. 
Consequently,  beach  states  do  not  change  through  the  full  range,  but 
rather  seem  to  alternate  between  the  longshore-bar-trough  state 
under  higher  breaker  conditions  and  the  rhythmic-bar-and-beach  state 
under  somewhat  reduced  breaker  heights.  Quite  simply,  it  appears 
that  breakers  do  not  reach  or  maintain  heights  large  enough  to  cause 
the  beach  and  surf  zone  to  become  fully  dissipative  and  they  do  not 
remain  small  long  enough  to  enable  the  crescentic  bar  to  weld  to  the 
beach  face. 

Despite  the  fact  that  the  morphodynamic  conditions  tend  to 
remain  arrested  between  the  longshore-bar-trough  and  rhythmic-bar- 
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and-beach  states,  these  states  are  among  the  most  mobile  (Short,  1981; 
Wright  1981)  with  the  result  that  profile  changes  are  appreciable  and 
rapid.  The  most  striking  changes  are  associated  with  the  onshore- 
offshore  migrations  of  the  bar.  Figure  47  shows  a  series  of  profiles 
measured  over  the  period  5-16  May,  1981.  During  an  earlier  phase 
of  this  period  (5-8  Kay)  comparatively  low  breakers  (H^Hm)  resulted 
in  onshore  bar  migration  at  a  rate  of  nearly  10  metres  per  day. 

However,  after  8  May,  onshore  migration  slowed  and  halted  owing  to 
a  return  to  higher  enegy.  The  associated  changes  as  they  appear  in 
plan  view  are  shown  in  Figure  25  from  which  it  can  be  seen  that 
longshore  migration  of  the  rip  and  of  crescentic  bar  ’’horns’’  accompanied 
onshore  bar  migration. 

From  Figures  25  and  47  it  is  evident  that  the  greatest  changes 
were  associated  with  the  bar.  The  beach  face  showed  much  smaller 
changes.  Most  of  the  beach  face  changes  which  occurred  on  a  given 
profile  were  actually  associated  with  the  alonghsore  migration  of 
rhythmic  features  rather  than  with  net  onshore-offshore  transfer  of 
material.  The  lower  mobility  of  the  beach  face  is  probably  due  to 
the  “filtering"  effects  of  the  bar  and  trough  system  in  maintaining 
a  relatively  constant  shore  break  amplitude. 

An  Atlantia  Coast  Analog 

Although  the  aerial  photos,  aerial  reconnaissances,  and  surveys 
of  Eastern  Beach  ocver  a  wide  enough  time  frame  for  us  to  make 
qualitative  conclusions  about  the  constancy  of  morphodynamic  state 
on  Eastern  Beach,  we  do  not  have  a  long  enough  time  series  of 
consecutive  surveys  to  ascertain  the  full  extent  of  baach  mobility. 

In  this  regard,  it  is  informative  to  make  use  of  results  from  an 
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analogous  example  in  a  similar  coastal  environment  elsewhere.  The 
analog  comes  from  the  vicinity  of  the  U.S.  Army  Coastal  Engineering 
Research  Center's  Field  Research  Facility  at  Duck,  North  Carolina 
on  the  Atlantic  coast  of  North  America.  Using  their  CRAB  surf-zone 
survey  vehicle,  the  Corps  of  Engineers  carries  out  frequent  routine 
surveys  of  the  beach  and  surf zone.  In  October-Nov ember ,  1981  two 
of  us  (L  D  Wright  and  P  Nielsen)  were  involved  in  an  experiment  at 
Duck  over  a  period  which  encompassed  a  series  of  storms  and  were 
impressed  by  the  morphodynamic  similarity  to  the  Eastern  Beach  site. 
Like  the  Gippsland  coast,  the  North  Carolina  coast  experiences 
frequent  storm  waves  which  typically  arrive  at  oblique  angles. 
Available  data  together  with  our  own  observations  suggest  that  the 
modal  morphodynamic  states  at  Duck  are  long shore-bar- trough  or 
rhythmic-bar- and-beach . 

Figure  48  shows  the  sweeps  of  superimposed  profiles  on  two 
transects  at  Duck  encompassing  a  period  of  11  months  and  including 
both  high  and  low  energy  extremes.  The  mobility  of  the  transects 
is  seen  to  be  confined  to  the  region  above  h  =  8  metres  and  is 
greatest  in  the  region  occupied  by  the  bar;  mobility  of  the  beach 
face  is  low.  The  bar  is  subject  to  rapid  migrations  over  a  horizontal 
distance  of  200  metres;  it  moves  offshore  during  storms  and  onshore 
under  diminishing  wave  conditions.  In  contrast,  the  beach  face 
apparently  advances  and  retreats  across  a  zone  only  ^20  metres  wide. 
These  characteristics  are  especially  evident  from  Figure  49  which 
shows  a  typical  low  energy  profile  and  a  post  storm  profile. 

Although  bar  positions  and  bar  crest  depths  are  vastly  different, 
the  beach  face  shows  only  modest  change  in  gradient  and  position. 
Furthermore,  it  is  readily  apparent  that  the  "cut"  of  the  beach  face 


Figure  49.  Comparison  of  a  typical  low  energy  or  "swell"  profile 
with  a  post-storm  profile  from  Duck,  N.C.  1981. 
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is  confined  to  the  subaerial  portion  above  a  shallow  step  and 
is  accompanied  by  accretion  on  the  segment  immediately  below  the 
step.  There  does  not  appear  to  be  significant  exchange  of  sediment 
across  the  trough  between  the  beach  face  and  the  bar.  An  important 
implication  is  that  not  only  are  the  bar  and  beach  face  regions 
dynamically  discrete  in  terms  of  their  dissipative  versus  reflective 
characteristics  but  are  also  at  least  quasi-discrete  in  terms  of 
their  mobility.  The  maintenance  of  standing  wave  nodes  in  the  trough 
may  provide  a  barrier  to  sediment  exchange. 

At  Duck,  the  sediment  population  was  strongly  bimodal  with 
medium  sand  composing  the  bar  and  coarse  to  very  coarse  sand  composing 
the  beach;  this  undoubtedly  influenced  the  modal  beach  state.  Similar, 
though  less  dramatic,  sediment  properties  prevailed  at  Eastern  Beach. 
Although  a  bimodal  sediment  distribution  is  not  a  prerequisite  for 
longshore-bar-trough  and  rhythmic  bar-and-beach  topographies,  it 
certainly  seems  to  favour  the  occurrence  of  those  states. 
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CONCLUSIONS 

In  approaching  coastal  problems  frcm  a  morphodynamic  point  of 
view,  we  are  concerned  with  the  complete  assemblages  of  depositional 
forms,  hydrodynamic  processes,  and  processes  of  sediment  redistribution; 
with  the  patterns,  sequences  and  extents  of  temporal  change  of  those 
assemblages;  and  with  the  relationships  of  the  assemblages  and  their 
temporal  variability  to  environmental  conditions.  In  the  foregoing, 
we  have  examined  a  Coastal  Boundary  Layer  system  in  the  presence  of 
environmental  conditions  characterize,  by  a  temporally-variable  deep¬ 
water  climate  subject  to  the  frequent  occurrence  of  obliquely-in- 
cident  storm  waves ,  shore  parallel  tidal  and  wind  driven  inner  shelf 
currents,  and  a  wide,  low  gradient  inner  shelf.  This  combination  of 
environmental  conditions  results  in  the  maintenance  of  beach  and  surf- 
zone  morphodynamic  states  characterized  by  longshore  bar-trough  and 
rhythmic  bar  and  beach  topographies  with  their  attendant  complex  suites 
of  fluid  motions.  We  have  attenuated  to  illuminate  some  of  the  con¬ 
nectivities  between  different  components  of  the  system  and  to  explain 
some  of  the  physical  mechanisms  which  cause  the  system  to  behave  as  it 
does. 

Specific  conclusions  and  implications  have  already  been  pointed 
out  in  the  appropriate  sections  and  subsections  of  this  report  and  do 
not  need  to  be  repeated  here.  However,  it  is  worth  highlighting  some 
general  points  of  possible  global  relevance: 

1.  Seaward  of  the  surf  zone,  tidal  currents  may  make  sub¬ 
stantial  secondary  contributions  to  the  shore  parallel  trans¬ 
port  of  sediment  even  in  microtidal  environments. 

2.  Over  rippled  beds  composed  of  coarse  material,  frictional 


dissipation  is  probably  appreciably  greater  than  has  been  con- 


ventionally  assumed.  Consequently,  the  configuration  and  roughness 
of  the  nearshore  profile  play  very  fundamental  roles  in  determining 
breaker  hight,  the  temporal  variability  of  breaker  heights,  and 
hence  the  modal  state  and  mobility  of  the  beach. 

3.  Although  longshore  bar-trough  and  rhythmic  bar  and  beach 
topographies  have  complex  process  suites,  incident  waves  remain 
dominant  right  up  to  the  beach  face  and  play  the  major  role  in 
suspending  sediments.  Incident  waves  also  exert  a  pronounced  in¬ 
fluence  on  the  structure  of  longshore  currents  by  contributing 
most  of  the  eddy  viscosity  near  the  bed.  Standing  waves  and 
possible  edge  waves  make  important  secondary  contributions  to  the 
net  drift.  Their  shore  normal  modes  may  locate  and  even  arrest 
the  trough  position  and  may  probide  a  weak  barrier  to  the  exchange 
of  sediments  between  the  bar  and  the  beach.  Cusps  and  rhythmic 
surf  zone  topography  are  probably  expressions  of  edge  waves  at 
different  frequencies. 

4.  The  bars  associated  with  longshore-bar- trough  and  rhythmic 
bar-and-beach  states  exhibit  high  mobility  characterized  by  rapid 
and  extreme  onshore-offshore  migrations.  However,  beach-face 
mobility  is  surprisingly  low  and  this  may  reflect  the  ability  of 
bar-trough  readjustments  to  partially  "absorb"  changes  in  breaker 
height  and  thereby  reduce  variations  in  shore-break  amplitude. 
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Appendix  A-1 :  Survey  Site  number,  location  and  shoreline  distance 
from  Cape  Howe  (N.S.W.  -  Victorian  border). 


EASTERN  VICTORIA 


No. 

Location 

Kms 

1 

Bastion  Point 

21 

2 

Mallacotta  (Golf  Course) 

24 

3 

Mallacotta  (Airport) 

27 

4 

Shipwreck  Creek 

32 

5 

Double  Creek 

35 

6 

West  Wingan 

56 

7 

Thurra  River 

77 

8 

Clinton  Rocks 

90 

9 

Bernm  River 

114 

10 

North  Conran 

130 

11 

Point  Ricardo 

142 

12 

North  Mario 

144 

13 

Corringal  Beach 

157 

14 

Ewings  Marsh  (Beach  Road 

165 

15 

Ewings  Marsh  (Pettman  Road) 

182 

16 

Eastern  Beach 

197 

17 

Barrier  Landing 

204 

18 

Bunga  Arm 

215 

19 

Ocean  Grange 

222 

20 

Stockyard  East 

235 

21 

Stockyard  Hill 

239 

22 

East  Causeway 

246 

23 

"Minnows" 

256 

24 

Paradise  Beach 

261 

25 

Golden  Beach 

266 

26 

South  Golden  Beach 

271 

27 

Carr  Creek 

276 

28 

Kullingal 

281 

29 

North  Seaspray 

286 

30 

Seaspray 

292 

31 

Denison  Creek 

293 

32 

McGaurans  Beach 

302 

33 

Woodside  Beach 

318 

34 

Reeves  Beach 

323 

35 

South  Reeves  Beach 

328 

Appendix  A-2 »  Frequency  of  occurrence  of  wind  speed  and  direction  at  Gabo  Island 
Lakes  Entrance  and  Tidal  River  (Source:  Bureau  of  Mecerology, 
Melbourne) . 
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.  >1  Hindi x  A-4 :  Direction  and  state  of  swell  and  sea  at  Gabo  Island  {Source:  bureau  of  Meteroloqy 
>!elbourne).  G*80  ISLAND 

DIRECTION  OF  SNELL  (1973-1982) 
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APPENDIX  B 

ESTIMATION  OF  FRICTION  FACTORS  AND  ROUGHNESS 
PETER  NIELSEN 


In  order  to  predict  frictional  dissipation,  we  must  be  able  to 

estimate  the  friction  factor  f  or  f  .  In  the  following  we  use  Jons- 

w  e  3 

son's  (1966)  definition  for  the  wave  friction  factor  f 

w 


T 

max 


1/20 fv"b 


,max 


(Bl) 


where  X  is  the  bed  shear  stress.  We  also  define  an  energy  dissipation 

factor  f  by 
e 

-  2  o 

tu,  =  —  pf  ir  (B2 ) 

d  3ir  e  D,max 

where  Xi^  is  the  time  averaged  energy  dissipation.  These  two  defini¬ 
tions  lead  to  identical  values  for  f  and  f  if  the  near  bed  velocitv 

e  w 

u^  is  simply  sinusoidal,  that  is  if 


^jmax 


sinuit 


and  x(t)  varies  as 


(B3) 


x(t)  =  x  Isinwtl  sinwt  (B4) 

max 

Lofquist  (1980)  has  measured  the  variation  of  x(t)  over  natural  wave 
ripples  and  found  that  it  can  differ  considerably  from  | sinwt |  sinwt, 
as  shown  in  Figure  Bl. 

Lofquist  also  derived  values  for  f  and  f  in  accordance  with 

we 


# 
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definitions  (Bl)  and  (B2) .  Figure  B2  shows  that  f  and  are  gene¬ 
rally  nearly  equal.  The  two  extreme  deviations  correspond  to  tests 
where  t ( t )  has  a  narrow  peak.  The  friction  factor  is  a  function  of 
the  Reynolds  number  t^u^  max,/,v  of  t*ie  roughness  to  semi  excursion 

ratio  X/a  (Jonsson  1966) . 
s 

When  the  flow  is  rough  turbulent  as  is  normally  the  case  under 
field  conditions  f^  (and  fg)  depend  only  on  a^/X. 

Several  authors  (Jonsson,  1966;  Kajiura,  1968;  Kamphuis,  1975; 
and  Grant  and  Madsen,  1982)  have  suggested  different  forms  of  the  re¬ 
lation  between  f^  (f^)  and  CL^/X  derived  in  analogy  with  steady  boun¬ 
dary  layer  flow.  According  to  Grant  and  Madsen  (1982)  this  analogy 

seems  reasonably  well  satisfied  when  a  /X  is  large  (a  /X  >  28)  as  in 

s  s  — 

Jonsson  and  Carlsen's  experiments,  however  the  variation  of  T (t)  as 
measured  by  Lofquist  shows  that  the  structure  of  the  boundary  layer 
over  large  natural  ripples  may  be  very  different  from  that  of  a  steady 
boundary  layer.  Figure  Bl  shows  that  t(t)  has  several  narrow  peaks  and 
troughs  that  are  related  to  the  formation  and  release  of  strong  ripple 
vortices.  A  boundary  layer  dominated  by  such  vortices,  generated  right 
at  the  bed  is  not  likely  to  be  logarithmic.  Therefore  it  seems  not 
worthwhile  to  compare  the  theoretical  merits  of  the  different  formulae 
in  great  detail.  Instead  we  will  use  the  most  handy  one  which  turns 
out  to  be  Jonsson's  in  the  explicit  form  given  by  Swart  (1974); 


fe  =  exp[5.213(^-)°-B4  -  5.977],  for  j-  <  .63 

s  s 


f 

e 


(B5) 


0.30 


(B6) 
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Figure  B3  shows  the  variation  of  f  with  A/a  . 

@  s 

The  upper  limit  of  0.30  is  rather  arbitrary  (Grant  and  Madsen 
argue  that  it  should  be  only  0.23) ,  and  it  is  traditionally  based 
on  Bagnolds  observation:  that  f£  did  not  increase  above  0.24  over  his 
artificial  ripples  in  sediment  free  water.  However, the  more  recent 
measurements  by  Carstens  et.  al.  (1969)  and  Lofquist  (1980)  using 
natural  sand  beds  indicate  that  f  can  be  0.30  and  even  larger,  and  it 
is  not  unlikely  that  the  presence  of  loose  sand  may  lead  to  friction 
factors  larger  than  the  maximum  values  for  fixed  beds.  This  is  suppor¬ 
ted  by  the  observation  of  Nakato  et.  al.  (1977)  that  the  presence  of 
loose  sand  enhanced  the  turbulence  intensity  by  a  factor  2  over  ripples, 
all  other  things  being  equal. 

The  influence  of  moving  sand  on  the  roughness  is  shown  in  Figure 
B4  where  we  have  plotted  measured  f^  values  from  Kamphuis  (1975)  and 
fg  values  from  Carstens  et.  al.  (1969).  Kamphuis  measured  the  shear 
force  on  a  flat  bed  of  fixed  sand  grains  and  Carstens  et.  al.  measured 
the  energy  dissipation  over  a  flat  bed  of  loose  sand.  We  see  that  the 
friction  is  nearly  one  order  of  magnitude  stronger  over  the  loose  sand 
for  the  same  values  of  a^/ Dg.  Dg  is  the  mean  grain  diameter. 

Kamphuis '  fixed-bed-measurements  are  well  predicted  by  equation 
(B5)  with 


A  =  2.5  D 

s 


(B7) 


For  the  loose-sand-data  the  corresponding  value  is  A  -  120Ds> 

When  the  bed  is  covered  by  ripples  or  other  bed  forms ,  these 
will  add  to  the  friction.  Following  Grant  and  Madsen  (1982)  we  will 
assume  that  the  combined  roughness  from  bed  forms  and  moving  sand  can 
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B4 :  Friction  factors  measured  over  flat  beds  of  fixed 

and  loose  sand.  We  see  that  the  friction  is  nearly 

an  order  of  magnitude  stronger  over  loose  sand  than 

over  fixed  sand  for  the  same  values  of  a  /D  .  The 

s  s 

curves  correspond  to  equation  B5  with  the  shown 
rouabness  lenaths. 


be  described  by 


A  =  (s  +  CM)F(e')Ds  +  K  ~  nr  (B8) 

where  is  the  ripple  height  ,  the  ripple  length,  9'  is  the  skin 
friction  Shields  parameter,  s  is  the  relative  density  of  the  sand 
(p  /p)  and  C  is  the  added  mass  coefficient  for  sand  grains  (=  0.5). 
Grant  and  Madsen  assume  the  value  28  for  the  constant  K  and  derive 

F(9  1 )  =  160  9c  [/S7  -  0.7/P]2  (B9) 

in  analogy  with  Owen  (1964)  who  studied  sand  grain  saltation  in  air. 
The  weak  point  in  this  analogy  is  that  the  vertical  length  scale  of 
the  sediment  motion  is  derived  from  the  idea  that  a  sand  grain  which 
hits  the  bed  with  horizontal  velocity  v.  may  bounce  off  vertically 
and  reach  a  height  of  the  order  of  magnitude  v?/2g.  This  is  probably 
true  in  air  where  the  relative  density  (sand  to  air)  is  of  the  order 
2000  and  inertia  and  gravity  may  dominate  over  the  drag  force,  but  in 
water  where  the  relative  particle  density  is  a  thousand  times  smaller 
a  sand  particle  will  only  move  a  distance  comparable  to  its  own  diame¬ 
ter  before  an  initial  upward  velocity  is  annihilatedby  the  drag  force. 
Furthermore,  sand  grains  do  not  bounce  off  in  water  because  of  the 
water's  cushioning  effect. 

However,  the  general  form  of  (B8)  is  physically  sound,  so  in  the 
following  we  will  derive  an  empirical  formula  based  on  (B5)  and  (B8) . 
The  skin  friction  shields  parameter  9'  is  calculated  from 
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T  f  UT 

w  D,max 

(s-1)  g  Dg 


(BIO) 


where  f  is  derived  from  (B5)  with  X  =  2.5D  .  First  we  look  for 
w  s 

the  value  of  the  constant  K.  Figure  B5  shows  observed  (via  equa¬ 
tion  (B5)) values  of  X/Cn^/L^)  versus  9'.  The  line  labelled  "Bag- 
nold"  is  derived  from  Bagnolds  (1946)  result 


L 

f  =  0.072  (^)  °-75 
s 


(Bll) 


which  can  be  transformed  into 


fe  =  3.75  (nJ/L .a.)0*75 


(B12) 


if  we  remember  that  the  ripple  steepness  in  Bagnold's  experiments  was 

n  /L  =  0.149,  and  that  f  =  f  /3  by  definition.  If  we  force  (B12) 

*  r  B  6 

to  coincide  with 


f  =  3.92  (-L)0.75  {B13) 

e  a 

s 

which  is  a  good  approximation  to  (B5)  for  0.1  <  —  <  0.5  we  find  that 

s 

Bagnold's  result  corresponds  to 


X  =  20.3  n  (B14) 

ii  r 


The  lines  labelled  "Jonsson  I"  and  "Jonsson  II"  refer  to  the  two  experi¬ 
ments  from  Jonsson  and  Carlsen  (1976)  where  XLr/n£  was  found  to  be  10.9 
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and  7.4.  As  8'  approaches  unity  XL ^/n2  tends  to  infinity  because 
nr/Lr  goes  to  zero  (e.g.  Nielsen  1981) .  The  optimal  value  of  K 
for  natural  sand  ripples  turns  out  to  be  approximately  8. 

Figure  B6  shows  (XQbg  -  8n£/Lr>/Ds  Plotted  against  6*.  Again 
the  scatter  is  very  considerable  but  there  seems  to  be  no  systematic 
deviations  between  different  grain  sizes.  The  flat-bed-data  from 
Carstens  et.  al.  conform  well  with  the  rest  although  they  correspond 
to  an  unstable  situation,  from  which  ripples  would  have  developed 
sooner  or  later.  We  see  that  the  formula 

(s  +  C  ) F (8 ' )  =  190/8 •  -  .05'  (B15) 

M 

gives  a  fair  representation  of  the  trend  while 

(s  +  C  ) F (8 1 )  =  504[/8T  -  0.16]2  (B16) 

M 

which  is  the  Grant  and  Madsen's  formula  with  C.,  =  0.5,  s  =  2.65  and 

M 

8c  =  0.05  seems  too  steep. 

We  will  therefore  recommend  the  following  formula  for  the  com¬ 
bined  roughness 

X  =  8n2/L  +  190/8'  -  0.05lD  (Bl7) 

r  r  s 

for  quartz  sand  with  s  z  2.65,  and  C,  z  0.5. 

M 

Figure  B7  shows  a  comparison  between  X/ag  derived  from  measured 
ffi  values  by  equation  (B5)  ,  and  values  predicted  by  equation  (B17) .  The 
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agreement  between  measured  and  predicted  values  is  acceptable  and 

the  agreement  between  the  corresponding  f  values  will  be  closer  be- 

,  d  f  ® 

1  e 

cause  — —  —  r— ■ — r-  is  considerably  smaller  than  unity  over  most  of 
f  d ( A/a  ; 

e  s 

the  practical  range. 

The  only  systematic  trend  seems  to  be  that  Carstens  et.  al.' s 
measurements  over  non  equilibrium  flat  beds  given  slightly  larger  f 
values  than  measurements  over  stabilized  rippled  beds  at  the  same  0' 
values.  This  may  mean  that  non  equilibrium  beds  are  rougher  because 
they  are  more  active  at  the  same  values  of  9'.  Savage  (1953)  obser¬ 
ved  similarly  that  the  friction  factors  were  larger  when  the  ripple 
pattern  was  not  in  equilibrium  with  the  wave  conditions.  Under  irre¬ 
gular  field  waves,  the  bedforms  can  never  be  in  equilibrium  with  the 
instantaneous  wave  conditions,  so  we  may  expect  slightly  larger  dis¬ 
sipation  factors  in  the  field  than  in  laboratory  experiments  with 
regular  waves. 

Figure  B8  shows  the  roughness  derived  from  Carstens  et.  al.'s 
"flat-bed-data"  versus  the  skin  friction  Shields  parameter.  We  see 
that  the  best  fit  for  these  data  alone  is 

(s  +  CJ  F (0 ' )  =  230/0'  -  0.05'  (B18) 

which  is  approximately  twenty  percent  larger  than  the  overall  best  fit 
given  by  equation  (B15) . 

In  conclusion  we  recommend  determination  of  f  (f  )  for  rough 

e  w 

turbulent  flow  over  natural  sand  beds  by  equation  (B5)  when  A/ag  < 

0.63  and  equation  (B6)  when  X/ag  >  0.63.  The  roughness  length,  X,  is 
best  estimated  by  equation  (B17) . 

For  prediction  of  the  ripple  geometry  under  field  conditions,  one 


J 


200 


230/e'-0-05 


Figure  B8:  Roughness  values  derived  by  equation  (B5) 
from  Carstens  et.  al,‘ s  "flat-bed-data". 
Note  that  only  four  highest  9' -values  cor¬ 
respond  to  fully  developed  rough  turbulent 
flow  according  to  the  criteria  given  by 
Jonsson  (1980). 
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may  use 


nr/Lr  =  °-342  ~  °-34  ^  for  6'  <  1. 


(B19) 


(B20) 


as  suggested  by  Nielsen  (1981).  When  9‘  >  1,  the  ripple  height  should 
be  taken  as  zero. 
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Non  dimensional  constant,  eg  (B8) 
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Relative  density  of  sand 
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Appendix  C-l:  continued 
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double  (inner-outer)  bar  systems.  6-13.12.8)  210-330ka  ground  survey 


/..tjundix  C-3 :  Selected  beach-surf  zone  morphological  and  sediment  characteristics  at  survey  sites 

LOCATION  BEACH  FACE  SURFZONE  TROUGH  TROUGH  X  CARB  BEACH  SURFZONE  OUTER 

GRADIENT  GRADIENT  WIDTH  DEPTH  0  BF  $  0  H.  H.  o.  M, 


Percent  Carbonate  Material 


Appendix  C-4:  Vegetation  tynes  at  survey  locations 
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depths  of  20m  and  10m  over  the  nearshore  profile  and  from  different  regions 
of  the  surf  zone. 

Over  the  nearshore  profile,  wave-induced  oscillatory  flows  dominate  to 
depths  of  20m  and  are  responsible  for  agitating  and  suspending  sediments. 
Despite  the  micro tidal  range,  reversing  shore-parallel  currents  and  super¬ 
imposed  wind  driven  currents  are  second  in  importance  and  determine  the  fate 
of  sediments  suspended  by  waves.  Frictional  dissipation  of  the  obliquely 
incident  waves  over  the  rough,  rippled  and  low-gradient  nearshore  profile 
is  particularly  pronounced  and,  as  a  result,  bottom  orbital  velocities 
measured  at  10m  depth  were  slightly  weaker  than  those  measured  at  20m.  Of 
fundamental  importance  to  surf  zone  processes,  nearshore  energy  dissipation 
by  bed  friction  causes  not  only  a  reduction  in  average  breaker  height  but 
also  significantly  reduces  the  temporal  range  of  breaker  conditions. 

The  morphodynamic  state  of  the  surf  zone  and  beach  remains  intermediate 
between  the  reflective  and  dissipative  extremes  and  is  characterized  for  most 
of  the  time  by  pronounced  longshore  bar-trough  or  rhythmic  bar  and  beach 
topography.  A  relatively  deep  (h  -  2m)  trough  separates  a  dissipative  bar 
from  a  steep  reflective  beach  face.  In  addition  to  the  oscillations  at  in¬ 
cident  wave  frequency  which  are  dominant  throughout  the  surf  zone,  important 
I  roles  are  played  by  surf  zone  standing  waves  at  subharmonic  and  infragravity 
frequencies  and  by  longshore  currents  and  rips.  Vertical  profiles  of  long¬ 
shore  current  velocity  within  the  trough  revealed  velocities  as  high  as 
0.25m  sec“l  to  within  0.075m  above  the  bed  and  were  not  at  all  logarithmic. 

The  profiles  are  explained  in  terms  of  the  large  near-bottom  eddy  viscosity 
caused  by  wave  induced  vortices  over  a  rippled  bed.  Direct  measurement 
of  vertical  profiles  of  suspended  sediment  concentrations  showed  that,  des¬ 
pite  the  presence  of  longshore  currents,  sediment  suspension  can  be  well 
modelled  in  terms  of  entrainment  by  waves  alone. 

Although  temporal  variability  of  deepwater  wave  conditions  is  very  large, 
the  substantial  reductions  in  breaker  variability  caused  by  increased  near¬ 
shore  dissipation  of  larger  waves  also  reduced  the  range  of  morphodynamic 
variability  of  the  beach  and  surf  zone.  The  longshore-bar- trough  and 
rhythmic-bar-and-beach  states  represent  the  modal  states  and  other  states 
rarely  develop.  However,  profile  mobility  is  considerable.  Most  of  this 
mobility  is  expressed  by  onshore-offshore  migrations  of  the  bar  and,  to  a 
lesser  extent,  by  longshore  migrations  of  rhythmic  features.  The  mobility 
of  the  subaerial  beach  face  is  comparatively  low  and  this  is  attributable 
to  the  relative  constancy  of  the  shore-break  amplitude .  The  Eastern  Beach 
modal  is  compared  to  an  analog  at  Duck,  North  Carolina,  U.S.A. 
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